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ABSTRACT 

We present a photometric survey of the stehar halo of the nearest giant spiral galaxy, Andromeda 
(M31), using Suprime-Cam on the Subaru Telescope. A detailed analysis of VI color-magnitude di- 
agrams of the resolved stellar population is used to measure properties such as line-of-sight distance, 
surface brightness, metallicity, and age, and these are used to isolate and characterize different com- 
ponents of the M31 halo: (1) the giant southern stream (GSS), (2) several other substructures, and 
(3) the smooth halo. First, the GSS is characterized by a broad red giant branch (RGB) and a metal- 
rich/intermediate- age red clump (RC). The /-magnitude of the well-defined tip of the RGB suggests 
the distance to the observed GSS field of (m — M)o = 24.73 ±0.11 (883 ±45 kpc) at a projected radius 
of ~ 30 kpc from M31's center. The GSS shows a high metallicity peaked at [Fe/H] > — 0.5 with 
a mean (median) of —0.7 (—0.6), estimated via comparison with theoretical isochrones. Combined 
with the luminosity of the RC, we estimate the mean age of its stellar population to be ~ 8 Gyr. 
The mass of its progenitor galaxy is likely in the range of 10^ to lO^M©. Second, we study M31's 
halo substructure along the north-west /south-east minor axis out to R ^ 100 kpc and the south-west 
major axis region at ~ 60 kpc. We confirm two substructures in the south-east halo reported by 
Ibata et al. (2007) and discover two overdense substructures in the north-west halo. We investi- 
gate the properties of these four substructures as well as other structures including the western shelf 
and find that differences in stellar populations among these systems, thereby suggesting each has a 
different origin. Our statistical analysis implies that the M31 halo as a whole may contain at least 
16 substructures, each with a different origin, so its outer halo has experienced at least this many 
accretion events involving dwarf satellites with mass lO^-lO^M© since a redshift of z ^ 1. Third, we 
investigate the properties of an underlying, smooth and extended halo component out to R > 100 kpc. 
We find that the surface density of this smooth halo can be fitted to a Hernquist model of scale radius 
~ 17 kpc or a power-law profile with cx 7^-2 i7±o.i5 contrast to the relative smoothness of the halo 
density profile, its metallicity distribution appears to be spatially no n- uniform with non-monotonic 
variations with radius, suggesting that the halo population has not had sufficient time to dynamically 
homogenize the accreted populations. Further implications for the formation of the M31 halo are 
discussed. 

Subject headings: galaxies: individual (M31) — galaxies: halos — galaxies: structure 



1. INTRODUCTION 

The advent of large ground-based telescopes and the 
Hubble Space Telescope {HST) has made it possible to 
resolve faint individual stars in nearby galaxies, thereby 
enabling us to study galactic archaeology for such exter- 
nal galaxies in the same way that one studies our home 
galaxy, the Milky Way. In particular, our knowledge of 
the stellar halo of the Andromeda galaxy (M31) has dra- 
matically expanded over the last two decades. Indeed, 
M31 provides us with a global yet detailed view of a 
galaxy with a similar morphology to the Milky Way. 
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The most significant discovery in the M31 halo is prob- 
ably its complex substructures such a s the giant south- 
ern stream (GSS) in the outer halo (Ibata et alj l2Q01l : 
Ferguson et^ 120021). This Andromeda GSS has at- 
tracted particular attention, as it reveals an on-going 
hierarchical formation process of M31 and thus places 
invaluable constraints on galaxy formation models. The 
GSS c ontains a high concentration of metal-rich stars 
(jibata et al. 2001; Guhathakurta et al. 2006) and is lo- 
cated behind M31 with a significant degree of elonga- 
tion along the l ine of sight ( McConnachie et al.i i2003i) . 
iFerguson et al.l (|2002l ) discovered significant substruc- 
tures around M31 by mapping red giant branch (RGB) 
stars with the Isaac Newton Telescope. 

The GSS was once thought to be tidal debris dis- 
persed by interactions with M32 or NG C 205 as Judged 
from its project ed trajectory on the sky ([Ferguson et al.l 
I2OO2I : [Choi et al...2002>) . but subsequent studies have re- 
vealed that it is probably not related to these two satel- 
lites (e.g., Ibata et al. '2004'; "McConnachie et aU I20Q4I : 
Guhathakurta et al. 2006; Kalirai et al. 2006a). In par- 
ticular, spectroscopic observations of GSS stars with 
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Keck/DEIMOS revealed that they show a markedly dif- 
ferent radial motion from M32 and NGC 205. Recently, 
A^-body numerical simulations have been performed in 
view of these kinematic results and these have showed 
that the GSS's progenitor had a highly e ccentric orbit 

and a mass of about 10^ 10^ M© (Fo nt et al.ll2QQ6al : 

iFardal et ani2QQ6l . l2QQ7t ). Furthermore, deep photomet- 
ric observations with ACS reaching down to b e- 
low the main-sequence turn-off (jBrown et alJ l2QQ6al |bl) 
detected intermediate-age stars in a GSS field located 
at a projected distance of R = 20 kpc from the M31 
center as well as in the general halo field at — ^P^ 
along the southern minor axis (jBrown et al.ll2QQ3h . 

In spite of these detailed studies however, it remains 
unclear what the detailed properties of the GSS's pro- 
genitor galaxy are. Besides, we do not fully under- 
stand the connection between the GSS and the other 
substructures in the M31 halo such as t he NE- and 
W-sh e lves, although the l atest simulations (jFardal et al.l 
I2QQ8I : iMori fc RichI l2QQ8l ) imply a connection among 
them. To unveil the origin of the GSS, it is of partic- 
ular importance to carry out the detailed analysis of its 
stellar populations by assembling a statistically signif- 
icant set of observational data such as color, metallic- 
ity, and age distributions of the stars, in addi tion to 
the other halo fields fe.g.. iBellazzini et al.l[2QQ3h . Stel- 
lar halos generally have such low surface brightness (SB) 
that detecting them beyond the Milk y Way is a ma- 
jor challenge. For example, Morrison (1993) estimate 
the SB of the Galactic halo at the solar radius to be 
jiiv ~ 27.7 mag arcsec"^. The M31 halo was long be- 
lieved to be an outward extension of its bulge, since sur- 
prisingly metal-rich stars (([M/H]) ^ —0.6) were found 
in the inner halo of M31 (e.g.. Mould & Kristian 1986; 
Couture et all 119951 : iHohand et al. 1996; Reitzel et al. 
19981 ). Furthermore, the SB profile measured along the 
south-east minor axis from integrated light and resolved 
star counts f ollows a de Vaucouleurs pr ofile out to 
R = 20 kpc (|Pritchet fc van d en Bergh"1994'). quite un- 
like the power-law behavior deduced for the halo of the 
Milky Way. Both the de Vaucouleurs profile and high 
metallicity are suggestive of an active merger history for 
the halo (or bulge) component of M31. 

The existence of metal-rich populations in the M31 
halo has b een confirmed by several subsequent studies. 
iDurrell eTaLI (|2QQll ) found that a high-metallicity popu- 
lation with [M/H] —0.5 is dominated in a field located 
at = 20 kpc along the south-east minor axis, based on a 
wide-field mosaic CCD camera on CFHT. They also dis- 
covered that 30%-40% of the stars at this location belong 
to a metal-poor populati on. Furth ermore, in their later 
complementary paper (Du rrell et al. 2004), they showed 
that a similar high-metallicity population is distributed 
in a field at = 30 kpc along the south-east minor axis. 
Therefore, they concluded that the outer halo shows lit- 
tle or no radial metallicity gradient. 

HST has made great contributions to the explo- 
ration of the inner h alo populations of M31. Using 
HST /WFFC2 [Bellazzmi et al. (2003) analyzed the data 
for 16 fields at different distances from the M31 center, 
from ~ 4.5 to 35 kpc, down to limiting V and / mag- 
nitudes of ~ 27. From the metallicity distributions in 
each field, obtained by comparison of the RGB stars with 



Galactic globular cluster templates, they detected simi- 
lar abundance distributions peaked at [Fe/H] ~ —0.6 in 
all of the sampled fields. Although in some fields a very 
metal-rich ([Fe/H] > —0.2) component is clearly present, 
possibly attributable to disrupted disk populations, they 
identified n o evident meta l licity grad ient in their lim- 
ited re gion. [Ferguson et al.l (|2005f ) and iRichardson et al.l 
(|2Q08l) investigated M31 inner substructures discovered 
bv lFerguson et all (|2002[ ). based on the deep HST /ACS 
survey spanning the range from R ^ 11.5 to 45 kpc reach- 
ing a few magnitudes below the RC. Their systematic 
surveys confirmed that the inner halo of M31 is strongly 
polluted by metal-rich and intermediate-age populations 
disrupted from the M31 disk and the GSS by interactions 
between M31 and dwarf galaxies. 

Brown et al.'s studies based on ultra-deep HST/ ACS 
observations have unraveled star- for mation histories in 
selected loca ti ons of the M 31 halo ([Brown et al.l [2QQ3I , 
I200I . [2006a ibl. I2OO7I . l2008h . With a minimum of 32 or- 
bits per field, these observations have been sufficiently 
deep to resolve stars below the oldest main-sequence 
turn-off. Brown et al.l (|2003f ) first demonstrated that a 
minor axis field projected at = 11 kpc contains an in- 
disputable population in age (predominantl y 6-10 Gyr) 
and abundance ([Fe/H] > -0.5). Brown etHI (|2006br ) 
present results from inner halo fields they associate with 
the "tidal stream," "outer disk," and "spheroid". In all 
cases, the fields were shown to have experienced an ex- 
tended star- formation history, i n co ntrast, th ey qua nti- 
fied differences between them. iBrown et al) (120071) in- 
vestigated a field at 21 kpc and found evidence that its 
population is marginally older and more metal poor tha n 
the inner halo field. Furthermore, iBrown et al.l (|2008l ), 
based on an additional field at = 35 kpc, found that 
the mean ages (the mean [Fe/H] values) at 11, 21 and 35 
kpc gradually change like 9.7 (-0.65), 11.0 (-0.87) and 
10.5 (-0.98) Gyr (dex), respectively. 

The inclusion of kinematic information has been ex- 
tremely useful but has also added another dimen- 
sion of complexity to t he puzzle. Spectroscopic stud- 
ies of iGuhathakurta e t al. (2005) uncovered success- 
fully the presence of an extended halo component hav- 
ing a power-law SB profile of R~^-^ over the range 
R = 30-165 kpc resembling the Milky Way's halo, 
while the SB of M31's spheroid at < 30 kpc has 
been known to show an R^^^ law resembling a ga lac- 
tic bulge (Pritchet fc van den Bergh 1994; Durr ell et al.l 
[2QQ1 : llrwin et all 120051 : [Gilbert et alIT2QQ6l ). They also 
reported that M31's halo may extend to R ^ 165 kpc, 
suggesting that halo stars of M31 and the Milky Way oc- 
cupy a substantial volume fraction of our Local Group. 
However, it remains to be seen whether this fiat portion 
of the halo profile i s affected by subs tructure. 

In the meantime iReitzel fc Guhathakurtal (|2002l ) an- 
alyzed a sample of 29 stars in a field at = 19 kpc 
on the minor axis and found the mean metallicity to be 
in the range of [Fe/H] = —1.9 to —1.1. Despite still re- 
maining uncertainties of calibration and sample selection 
issues, their estimated mean metallicity is significantly 
lower than th e results dedu c ed frorn the above photomet- 
ric analyses. iKalirai et al.l ([2006bl) and I Chapman et al.l 
(2006) recently discovered a metal-poor halo out to = 
165 kpc that is similar to that of the Galaxy, both in 
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terms of its metallicity, radial profile, as well as the 
amount of substructure. 

As noted above, the stellar halo of M31 as well as the 
Milky Way is filled with quite complex fossil records. 
Our principal goal is thus to identify each component 
of the halo and clarify its formation history. To do so, 
we need to obtain both the global and local structures 
in the stellar halo and deduce plausible origins of them. 
Detailed studies based on large numbers of stars, such 
as with wide-field imagers, are required to investigate 
the fundamental nature of the large halo regions in M31. 
Indeed, our external viewpoint for this gala xy makes it 
possi ble to study the entire area of its halo (jlbata et al.l 
|2QQ7[ ): detailed stellar maps will allow us to detect giant 
streams as well as a wealth of fainter substructures at 
SB levels of ~ 30 mag arcsec"^. In the south-quadrant 
halo of M31, Ibata et al. also discovered much fainter 
stream-like substructures reaching ~ 34 mag arcsec"^ as 
well as previously unknown surviving dwarf satellites and 
peculiarly extending outer globular clusters. The prop- 
erties of these faint substructures in the outer halo, such 
as radial distributions of surface brightness and metal- 
licity, have been reproduc ed by s ome recen t numerical 
simulations (e.g., i Bullock fc Johnstonll2QQ5l : lAbadi et al.l 
I2QQ6I : iJohnston et al. 2008), which thus provide us with 
useful theoretical templates to understand the formation 
history of the stellar halo. 

Here we report on our wide-field, photometric obser- 
vations of the M31 halo, using Subaru/Suprime-Cam, 
aimed at obtaining the detailed physical properties of the 
M31 halo, such as color-magnitude diagrams (CMDs), 
metallicity distributions, surface brightness and age, 
comp ared with previou s studies using CFHT/MegaCam 
(e.g., 'Ibata e t aJ] 120071) . Our Suprime-Cam survey is 
deeper and covers as-yet unexplored fields along the 
north-west minor axis of the halo as well as some fields 
in the south-quadrant halo already explored by other 
telescopes. A large set of our imaging data provides an 
important clue to understanding the complex halo sub- 
structures and their origin in the course of the formation 
history of the M31 halo. 

The layout of this paper is as follows. In § [21 we 
present observations and detailed procedure for calibrat- 
ing the photometric Suprime-Cam data. § [3] is devoted 
to our results for the quantitative analysis of the GSS, 
including the morphology of the obtained CMDs, the dis- 
tance estimate using the tip of the RGB stars (TRGB), 
their metallicity distribution using the observable globu- 
lar cluster templates and theoretical isochrones, and their 
age distribution based on the luminosity distribution of 
the stars. In § IH we expand our method to the entire re- 
gions of the M31 halo in order to globally determine the 
fundamental properties of the M31 halo. Then, we elab- 
orate the technique for separating the foreground and 
background contaminants. The spatial distributions of 
the stellar populations in the M31 halo are also presented 
in the section. Furthermore, we present the CMDs of the 
detected streams and other spatial substructures, and 
show the radial profiles and the radial metallicity distri- 
butions of the stellar populations in the halo. Finally, we 
discuss the implications of our results and compare with 
previous studies. In §[8] we draw conclusions and future 
directions. 

Throughout this work we adopt the M31 Cepheid dis- 



tance mod ulus of (m - M)n 24 .43 ± 0.06 (- 770 kpc) 
from Free dman fc Madorq (|1990[ ) , if not otherwise spec- 
ified. Reddening corrections are appli ed in ea c h fiel d 
based on the extinction maps of Schle gel et aP (UoOSh . 
and the iDean et all ()1978l) reddening law E{V - I) = 
1ME{B - V) and Aj = 1.31£;(V - /). We also adopt 
the convention that R and r denote, respectively, the 
projected radius and three-dimensional radius from the 
M31 center. 

2. DATA 
2.1. Suprime-Cam Observations 

In our observational studies of Andr omeda's stellar 
halo, we use the Suprime-Cam imager (iMivazaki e t al.l 
I2002f ) on the 8.2-m Subaru Telescope (lye et al. 200^ 
on Mauna Kea in Hawaii. Suprime-Cam consists of ten 
2045 X 4070 CCDs with a scale of 0'.'202 per pixel and 
covers a total field-of-view 34' x 27'. Using this wide-field 
imager, we have carried out a systematic imaging survey 
of two series of fields far from the M31 center: fields 
along the minor and major axis of the spheroid (here- 
after referred to as "halo" fields) and fields along the 
Andromeda GSS ("stream" fields). Our targeted fields 
of the halo and stream fields are, respectively, located at 
between 22 and 91 kpc (Figure [T]) and at - 32 kpc (2?4) 
from the M31 center. Tables [T] summarizes the names of 
these targeted fields and their observations. 

During four nights in August 2004, we obtained the 
Suprime-Cam image of the one stream field (hereafter 
referred to as GSSl) and also made imaging surveys for 
eleven halo fields in the south part of M31 (SEl-9 and 
SWl-2). In addition to these targeted fields, we obtained 
the imaging data at the control field (hereafter referred 
to as CompL103Bm21) to estimate the number of fore- 
ground dwarf stars of the Milky Way and unresolved 
background galaxies in each M31's field. The weather 
conditions in August 2004 were unfortunately slightly 
poor, with seeing varying from 0'.'7 — 1%. Therefore, 
we carried out additional imaging of GSSl, SEl, SE6, 
SE9 and SWl on a photometric night in August 2005 in 
order to calibrate the data. In August 2007 and August 
2008, the halo fields in the north part of M31 (NWl-15) 
were observed in photometric conditions, along with the 
control fields (CompL103Bm25 and CompL103Bml6_5). 

The observations were made with Johnson V-band and 
Cousins /-band filters. Typical exposure times of each 
object field observed in 2004-2005 and 2007 are 1800 sec 
and 720 sec in V and / band, respectively, and those in 
2008 are 1440 and 2160 sec in V and / band, respectively, 
so that the images reach down to magnitude of the RC 
of M31. Each field was observed with adequate dither- 
ing pattern to cover gaps between adjacent CCD chips. 
Some of these fields were observed in short exposures 
of 5 and 10 s (12 and 15 s) in V (Ic) bands. They are 
used to improve photometric accuracy and to investigate 
space variation of foreground dwarf stars along the galac- 
tic latitude at the same galactic longitude. In addition, 
all the targeted fields in the M31 halo are overlapping 
with neighboring fields in order to calibrate the adjacent 
fields whose photometric zero points are unknown in poor 
condition, such as SE2, SE3 and so on. 

The position angle of each field observed in 2004-2005 
is 0°, i.e.. North is up and East is to the left of each 
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image. In the 2007-2008 runs, the position angle of each 
field is 38°, i.e., the longer side of each image is aligned to 
the minor axis of M31 in order to obtain many pointings 
along the far outer part of M31's halo. For more details, 
see Tabled! 

2.2. Reduction and Photometry 

The raw data were reduced in the standard procedures, 
with the software package SDFRED, a useful pipeline 
developed to opt i mally deal with S uprime-Cam images 
(|Yagi et al.l[2QQ2l : IQuchi et al.l[200l . Using SDFRED, 
we obtained a mean-stacked image for each field, which 
covers 34' x 27' field-of-view. 

We then conducted PSF-fitting photometry u sing the 
IRAF version of the DAOPHOT-E software (Stetson 
[1987). We adopted a 3 cr detection threshold for the 
initial object detection/photometry pass and repeated 
the PSF-fitting photometry twice with 5 a and 7 a de- 
tection thresholds, respectively, in order to account for 
blended stars. A stellar PSF template was constructed 
from 10-100 bright, not saturated, isolated stellar objects 
per image. Finally, we merged two independent V-band 
and /-band catalogues into a combined catalogue using 
2-pixel matching radius. Note that we cannot detect ob- 
jects near the edge of a field in case of adopting a 1-pixel 
matching radius. 

When conducting photometry, we divide a Suprime- 
Cam field into a few subfields and photomery is carried 
out for each subfield separately. For example, for fields 
with PA= 0° we divide them to four subfields (labeled 
by "a", "b", "c" and "d"), while for fields with PA= 38°, 
we divide them to two subfields (one is not labeled and 
the other is labeled by "d"). If we perform detection and 
photometry using DAOPHOT without dividing to sub- 
fields, stellar-like objects near the edge of field-of-view 
(34' X 27') are not properly detected because the actual 
PSF is largely varying depending on the location. Even if 
we consider the space variation of PSF by means of mak- 
ing a model PSF function to detect objects {PSF task), 
it results in inducing erratic values of sharpness and 
chi in some DAOPHOT parameters. Actually, when we 
have detected objects without dividing one field-of-view, 
we have failed to detect approximately 20% of stellar 
samples of M31 red giants and Milky Way dwarf stars 
ident ified b y Keck/DEIMOS spectroscopic observations 
(see g223j). 

2.3. The Artificial Star Experiments 
2.3.1. Completeness 

To evaluate incompleteness due to low signal-to-noise 
ratio and crowding. We have conducted the artificial star 
experiments as detailed below. We have added a suffi- 
cient number of artificial stars to all the original images, 
where each star has the PSF as determined from many 
bright stars in each frame in § 12. 2[ and the range of the 
magnitudes is 20 < V, / < 28 with a binning step of 
0.2 mag. Note that for a portion of the large images the 
binning is 0.5 mag in the brighter parts with V = 23 mag. 
To prevent these stars from interfering with one another, 
we have divided the Suprime-Cam frames into grids of 
cells of 200 pixels width and have randomly added one ar- 
tificial star to each cell for each run. In addition, we con- 
strain each artificial star to have 20 pixels from the edges 



of the cell. We note that a similar procedure to our s has 
recently been adopted by iPiotto fc Zoccalil (|1999[ ) and 
iBellazzini et al.l (|2002[ ). In this way we can control the 
minimum distance between adjacent artificial stars. At 
each run the absolute position of the grid is randomly 
changed in a way that, after a large number of experi- 
ments, t he stars are unifo rmly distributed in coordinates 
(see also lTosi et all 1200 ID . 

The stars, which have a Penny 2 function expressed as 
the analytic component of the PSF model computed by 
the DAOPHOT-E PSF task and are quadratically vari- 
able over the image, were added on the original frame 
including Poisson photon noise. Each star has been 
added to the V and / coadded frames. For the artifi- 
cial frames we have performed detection and photome- 
try following the method described in ? l2.2l and extracted 
stellar-like objects by applying our established selection 
criteria based on DAOPHOT parameters (see ? I2.3.2|) . 

About 500 to 4300 stars, were configured for each run, 
and the total number is about 18,000 to 140,000 in each 
band image, where the number of added artificial stars 
depends on the subfield size. After the above-mentioned 
detection/photometry and selection, the final recovered 
objects were regarded as stars, if they reside within 2 
pixel of the position of the added star and their magni- 
tudes are within 0.75 mag of the assigned values. In this 
task, the completeness is defined as a number fraction of 
recovered stars in added stars. In Figure [21 we plot the 
completeness fractions (= A^output/Mnput) as a function 
of input V (filled circles) and / (open rectangles) mag- 
nitudes in the NW3 field as estimated by the current 
artificial star experiments. Note that the continuous line 
is not a fit to the data but is simply the line connecting 
them. Alphabetical letters labeled in each field of the 
south quadrant area of the M31 halo, such as "a", "b", 
"c" and "d" represent the subfields which are obtained by 
dividing a single Suprime-Cam field into four subfields. 
In these subfields, the 25% overlapping regions (labeled 
by "c" for SE and SW fields and by "d" for GSS fields) 
indicate slightly deep completeness limits. In contrast, 
the regions which are overlapping the nearby fields along 
the North- West minor axis of the M31 halo are labeled by 
"d" , meaning "deep" . The rightmost column of Table [T] 
shows the 50% completeness limits of V- and I-band in 
each subfield. In all the fields the completeness is larger 
than 80% at / < 24.5, and our imaging data are much 
deeper than t he previous studie s using grou nd-based tele- 
scopes., e.g;. iMcConnachie et a l. (2003), iDurrell et al.l 
(|2QQll . l2QQ4f ) andll bata et all (|2007 ). The 50% complete- 
ness is reached dit V > 25.70 and / > 24.85 in all the 
fields. 

2.3.2. Selection of Stellar-Like Objects Based on 
DAOPHOT Parameters 

In extracting stellar-like objects from the catalogs 
containing all the detected objects, we apply selection 
criteria based on DAOPHOT parameters such as the 
sharpness, chi, and the photometric error, which are 
simulated in the above artificial star experiments. The 
sharpness parameter estimates the intrinsic angular size 
of the measured objects. It is roughly defined as the 
difference between the square of the width of the object 
and the square of the width of PSF and has values close 
to zero for single stars, large positive values for blended 
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TABLE 1 

Observational status for target fields 



Field 




Date 


Filter 


J--J-A.L^UO Lll 
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Qppi n 

kJC^C^lll^ 






1 nr\r\r\ / nonr»r» 1 
\^2\J\J\J / ^ZUKJU J 


mm/ WW 
111111/ y y y y 




(sec) 




^^iiidgy 




GSSl 


00h47™38?ll 


08/2004 


V 


600 X 5 


1.7 


1('07 


25.95 


32 kpc) 


+39°05'57'/4 


08/2004 


I 


300 X 6 


1.2 


0('92 


25.00 


SEl 


00^50™09?16 


08/2004 


V 


600 X 3 


1.3 


1^00 


26.40 


24 kpc) 


+40°08'36'/4 


08/2004 


I 


300 X 6 


1.1 


0^71 


25.66 


SE2 


00^51"^38?06 


08/2004 


V 


600 X 5 


1.2 


1^13 


26.34 


29 kpc) 


+39°55'03'/4 


08/2004 


I 


300 X 6 


1.1 


0^93 


25.24 


SE3 


00h53"^06?96 


08/2004 


V 


600 X 3 


1.4 


1^03 


26.30 


34 kpc) 


+39°41'30'/4 


08/2004 


I 


300 X 6 


1.1 


0^82 


25.40 


SE4 


00h54"^35?86 


08/2004 


V 


600 X 5 


1.3 


l'/13 


26.34 


39 kpc) 


+39°2r5r/4 


08/2004 


I 


300 X 6 


1.1 


0'/87 


25.27 


SE5 


00h56"^04?76 


08/2004 


V 


600 X 5 


1.5 


V/21 


26.24 


44 kpc) 


+39° 14^24^4 
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Fig. 1. — Locations of our Subaru/Suprime-Cam fields (rectangular areas), overlaid on the color-composited stellar density map [taken 
from Figure 53 of the Ibata et al. (2007) paper], in which red, green, and blue colors show, respectively stars with —0.7 < [Fe/H] < 0.0, 
— 1.7 < [Fe/H] < —0.7 and —3.0 < [Fe/H] < —1.7. Our Suprime-Cam fields characterized by blue, yellow, and green rectangles were 
obtained in 2004, 2007 and 2008, respectively. Table 1 presents a log of these observations. 
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doubles and partially resolved galaxies, and large nega- 
tive values for cosmic rays and blemishes. The niter pa- 
rameter is the number of iterations the solution required 
to achieve convergence. The chi parameter assesses the 
estimated goodness-of-fit. It is the ratio of the observed 
pixel-to-pixel mean absolute deviation from the profile 
fit, to the value expected on the basis of the noise as de- 
termined from Poisson statistics and the readout noise. 

The simulation results of the four parameters to select 
stellar-like objects are illustrated in Figure [3l The solid 
and dashed lines show mean value and 4 a deviations of 
each parameter measured for artificial stars, respectively. 
We choose as "good" stars those sources lying within two 
hyperbolic envelopes around the sharpness value of zero 
(Fig. [3] bottom left panel), that is, within the 4 sigma 
deviation drawn by the dashed lines, below the 4 sigma 
lines of the niter value (top right panel) and the chi value 
(bottom right panel), and to the right of the maximum 
allowed a error at a given magnitude (top left panel). 
The remaining cosmic-ray blemishes and resolved back- 
ground galaxies were then rejected. However, this sigma 
criterion slightly changes with the sky condition because 
this simulation is prone to misclassify nearly point-like 
background galaxies as stellar objects under the worse 
seeing condition. We have tested the seeing effect by 
comparing the overlapping region between the NW7 field 
with good seeing and the NWS field with bad seeing. To 
adjust the surface brightness of the NWS subfield to that 
of the NW7 subfield, we have found a 1 sigma shift of 
our adopted selection criteria. Therefore, we apply the 
3 sigma criteria for the poor data with the larger seeing 
size than OVS. 

However, this morphological segregation method can- 
not completely distinguish M31 halo stars from other 
point sources such as Galactic foreground dwarf stars, 
and compact extragalactic objects. Thus, we have to re- 
move these contaminations from the targeted fields sta- 
tistically. To estimate the contaminations, we adopted a 
control field located at the Galactic latitude comparable 
to that of our targeted fields, on the assumption that it 
has almost the same abundance of foreground and back- 
ground objects as that in the M31 spheroid fields. Actu- 
ally, we found uniform galaxy distribution in our control 
fields. As described in § 14.1. If to further refine our con- 
trol field, we will consider the finite spatial gradient in 
the number of foreground disk dwarf stars of the Milky 
Way along the Galactic latitude. 



2.3.3. Photometric Errors 

These artificial star experiments are implemented to 
ascertain not only the completeness but also more ac- 
curate photometric errors in our data. These are quan- 
tified as the average of differences between input and 
output magnitudes of an artificial star at a specific mag- 
nitude, shown in Figure [H The mean magnitude differ- 
ence Amag = magij^py^ — mag^^^p^^ is consistent with 
for magnitudes brighter than a S0% completeness limit 
{V - 26 and / - 25.5 in case of the NW3 field). The 
shift toward systematically brighter measured magni- 
tudes, which signals the occurrence of blending of the 
stellar images, is smaller than 0.09 mag at the magni- 
tude bin corresponding to 50% completeness level. The 
right panels of Figure [1] show mean error as a function of 
magnitude as measured from these simulations. 

2.4. Comparison with Published Data 

2.4.1. Detection featuring Keck/DEIMOS Spectroscopic 

Samples 

In SE3 field, there are som e M31's RGB stars iden- 
tified bv [Gilbert eFall ([2QQ6h . Using the DEIMOS in- 
strument on the Keckn 10m telescope, their RGB sam- 
ples have been confirmed based on five criteria: radial 
velocity, photometry in the intermediate- width DD051 
band to measure the strength of the MgH/Mg6 absorp- 
tion features, strength of the Nal AS190 absorption line 
doublet, location within an (F — /, /) color-magnitude di- 
agram, and comparison of photometric (color-magnitude 
diagram based) versus spectroscopic (Call AS500 triplet 
based) metallicity estimates. Therefore, to test the va- 
lidity of our detection method, we match up our photo- 
metric catalogs with their spectroscopic catalogs using a 
10-pixel matching radius. 

Figure [5] shows a small section of the SE3 field (2250 x 
5200 pixels, the dotted rectangle in Figure [6|). This 
field is located close to the center of the DEIMOS 
spectroscopic field, at a2ooo = 00^54"^0S^34, ^2000 = 
+39°41'5lV7. The colored circles denote those objects 
whose spectra are measured by Keck: M31 red giants 
and Milky Way dwarf stars (red, 61 objects), background 
galaxies/QSOs (cyan, 31 objects), alignment stars (blue, 
6 objects) used to carry out fine astrometric alignment of 
the DEIMOS multislit mask, spectroscopic failures (yel- 
low, 5 objects) and guide stars (magenta, 2 objects) only 
used for coarse astrometric alignment - no spectra are 
obtained for them. In contrast, green X points show our 
Subaru objects corresponding to all the Keck samples in 
our final star catalog. 

The number of stellar-like objects detected by our Sub- 
aru observation is 56 objects, and the number of missing 
objects is less than 7% of the total. Note that one out 
of 61 spectroscopic stellar-like objects is in the masked 
region, and therefore the total number is 60 objects. 
Figure [3 shows the magnitude difference between these 
different systems. The 4 missing objects are removed 
from the final catalog because they have somewhat large 
sharpness and chi values in their magnitude. 

We have also detected 12 objects classified as back- 
ground galaxies/QSOs by the Keck spectroscopic ob- 
servation. This is 40% of the total number, and the 
contaminations attributable to these background objects 
are moderately removed. This suggests that all the 
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Fig. 3. — Photometry selection criteria based on the photometric error, niter, sharpness and chi parameters from DAOPHOT programs 
are plotted as functions of V and I magnitudes. All of the detected objects after matching up to separate V and I catalogs derived from 
the observation of NW3 field are illustrated as light gray dots. The solid lines show returned mean values of each simulated parameter, 
and the dashed lines denote the 4 sigma deviations. We select final objects within the 4 sigma deviations. 



background contaminations are never rejected using only 
DAOPHOT morphological parameters. Therefore, the 
remaining background contaminations are removed in 
comparison to the data of Control Fields located outside 
M31's halo (for more details, see § I2.3.2p . 

2.4.2. Photometry featuring KPNO Telescope Samples 

In this section, we compare our photometric data with 
other photometric data taken with the Mosaic camera on 
the Kitt Peak National Observatory (KPNO) 4m tele- 
scope in the Washington System M and T2. Pho t omet - 
ric transformation relations from^ Majewski et alJ (|2QQQl ) 
were used to derive Johnson- Cousins V and / magnitudes 
from the M and T2 magnitudes. The color transforma- 
tion follows: 

{V -I) = {V - T2) = -0.006 + 0.008(M - T2). (1) 



With the assumption I = T2 and the above relation, we 
converted to a F magnitude 

V = M - 0.006 - 0.200(M - T2). (2) 

Figure [6] shows sky positions observed by three differ- 
ent telescopes, which are used to evaluate the efficiency 
of our detection/photometry. The "aO" field is associ- 
ated with our SEl, SE2 and SE3 fields, and the "mask4" 
field covers our SE3, SE4 and SE5 fields, and the "m6" 
field overlaps our SE8 and SE9 fields. 

We compare the magnitudes based on the Subaru cal- 
ibration with the magnitudes based on the KPNO cali- 
brations after matching up two different catalogs using 
a 3-pixel matching radius. We extract stellar-like ob- 
jects based on the DAOPHOT parameters {merr < 0.02, 
—0.3 < sharp < 0.3 and chi < 4). In Figure [71 we 
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Fig. 4. — Typical mean error as a function of magnitude for NW3 field. On the left we plot difference between the measured and input 
magnitude as a function of input magnitude for all the stars from our completeness simulations, while the panels on the right show mean 
error as a function of magnitude as measured from these simulations. Only those simulated stars that were recovered within 2 pixels of the 
input position and that had magnitudes within 0.75 of the input magnitude are considered to be recovered. There is no systematic shift of 
the mean magnitude difference, and the typical 1 a errors at the 50% completeness magnitude limits are of the order of 0.3 mag. 



present a comparison of the two independently deter- 
mined V- and /- band flux measurements for the three 
KPNO fields. The fact that the number of faint stars 
with V ^ 23.5 and / ^ 22.5 decreases with decreas- 
ing fiux is due to low signal-to-noise in fainter part of 
the KPNO survey. The individual Fsubaru (^Subaru) and 
^KPNO (^KPNo) are found to be in good agreement with 
one another over most of the magnitude range, indicating 
that there are no systematic variations in our magnitude 
measurements. 

3. THE ANDROMEDA GIANT STREAM 

3.1. Morphology of the Color-Magnitude Diagram 

Figure El shows the CMDs for the GSS (left) and for 
the control field (right). In this section, we use the SE3 
field as the control field. The clearest feature in the 
stream field is a broad RGB characterized by low-mass, 
hydrogen-shell-burning stars. It also follows that the 
GSS has a broad metallicity distribution and to some ex- 
tent metal-rich R GB attributable to a lar ger stellar enve- 
lope opacity (e.g. JSalaris fc Cassisil[2QQ5h. in ps;ood agree- 
ment with the previou s studies (e.g.." McConnachie et al.l 
I2QQ3I : iFerguson etHI .2005. ). The CMD of the control 



field also shows an apparently broad RGB feature, al- 
though it is less dense than that of the GSS. It is worth 
noting that our control field is located slightly more re- 
m ote 5 kpc) from M 31's center than the field analyzed 
by I Purrell et al.l (12004]) on the minor axis, thereby sug- 
gesting that metal-rich halo stars are distributed even i n 
the outer parts of the M31 halo (e.g. Jlrwin et al.ll2005f ). 

To extract the detailed distribution of stellar popula- 
tions in the GSS alone, the field stars other than the 
stream stars must be removed from the CMDs. For 
this purpose, we have made statistical subtraction of 
the field stars taken from the CMD of t he co ntrol field, 
following the method bv iDurrell et all (|2004f ). In Fig- 
ure [9l we show a colored contour for the log-scaled CMD 
of GSSl field after this procedure, illustrating detailed 
characteristic features of the GSS stellar populations (cf. 
Bellazzini et al. 2003). It follows from this figure that in 
addition to the RGB feature, we identify a Red Clump 
(RC) at / ~ 24.5, structure in the CMD consisting of a 
number of helium core burning stars. More accurate es- 
timates of the RC magnitude will be made below, by cor- 
recting for the incompleteness of magnitudes (Sec. 13. 4|) . 
The detection of a clear RC thus indicates the presence 



10 



Tanaka et al. 




Fig. 5. — A small section of the SE3 field. This is a V-band 
coadded image with the total exposure time of 30 min and covers 
a field-of-view 7.6' x 17.5^ 



of metal-rich and/or young populations in the GSS. 

In Figure [H there is also a second peak at / ~ 23.5, 
which is apparently brighter than the RC. This fea- 
ture can be regarded as an Asymptotic Giant Branch 
(AGB) bump, not a RGB bump. This is for the rea- 
son that while most stars in the GSS appear to be old 
and/or metal-rich populations, a RGB bump with such 
old and/or metal-rich populations lies at a fainter magni- 
tude than a RC magnitude in the CMP , following various 
theoretical models as calibrated by Alves fc Sarajedinil 
(1999.) . There is also a possibility that a clumpy fea- 
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Fig. 6. — Sky positions observed by three different tele- 
scopes, which they are used to evaluate efficiency of our Detec- 
tion/Photometry. The gray rectangles represents the position and 
area of the Subaru/Suprime-Cam image. The solid squares indi- 
cate the approximate positions of the KPNO fields. The dotted 
vertically long rectangle approximates the size and position angle 
of the DEIMOS spectroscopic slitmask. 



ture being slightly brighter than a RC is just a portion 
of the red horiz o ntal b ranch (RHB) clump as shown in 
iBellazzini efall (j2QQ3f ). Their G76 field having such a 
feature is located at an edge of a luminous thin disk of 
M31, so only young, metal-rich populations like disk stars 
give rise to this kind of a feature slightly brighter than a 
RC magnitude. 

At / ~ 20.5, we also recognize the tip of the RGB 
(TRGB). The presence of a few stars brighter than the 
TRGB magnitude level presumably suggests the presence 
of long-period variable stars in metal-rich and/or young 
population like bri ght AGB stars, as is the case in the 
halo of NGC 5128 dReikuba et al.l[2003h . 

3.2. Distance to the Andromeda Giant Stream 

As mentioned above, the TRGB is a useful in dicator to 
estimate the distance to res olved galaxies (e.^.ALee et al.l 
[1993; Madore & Freedma3 119951 : iFerrarese etld] I2QQQI : 
ISalaris & Cassisii |2QQ5.) . If there are a sufficient num- 
ber of old and metal-poor RGB stars in a targeted 
field, the TRGB is easily detected as a sharp cut- 
off of the luminosity function (LF) with the applica- 
tion of an edge detection algorith m, Sobel filter with 
a kernel [-1, -2 0, -h2 , +1] (Mador e FreedmanI [19951 : 
ISalaris fc Cassisill2QQ5l ). In this manner, we have identi- 
fied the magnitudes of the TRGB from /-band LFs and 
determined the distance to the GSS. The TRGB is safely 
detectable, since it is bright enough not to be affected 
by incompleteness (i.e., more than 95% completeness for 
both V and /-band magnitudes near the TRGB) and 
since it is fainter than a saturation level (/ ~ 19.5) of 
the CCDs. 

Detection of the TRGBs for the GSSl field is shown 
in Figure [To] (see also Fig. [Ej). These smoothed LFs are 
constructed from the raw catalog by treating each star 
as a unit Gaussian with width <j(/) and summing all the 
Gaussians (Durrell et al...2QQl.) . Taking into considera- 
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Fig. 8. — Color-magnitude diagrams in the /q vs. {V — I)o plane for GSSl and the Control Field (SE3). The dashed lines denote the full 
ranges of the 50% completeness levels for the mosaic images of the Suprime-Cam. Extended background galaxies have been removed via 
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Fig. 9. — A log-scaled CMD for GSS 1. The peak number densities 
of the AGE bump and the RC amount to 224 and 657 stars per 0.07 
X 0.07 mag box, respectively. Stellar density increases with redder 
color. The main features are distinguished by an edge of stellar 
density: the TRGB at Iq ~ 20.5 and by the reddest peaks of high 
stellar density: the AGB bump at Iq ~ 23.5 and the RC at Ip ^ 
24.5. The solid lines are theoretical isochrones (VandenBerg ^ al.l 
l2006t ) of age 12 Gyr spanning the metallicity range [Fe/H]= —2.14, 
-1.41, -1.14, -0.71, -0.20 and 0.00. The dotted line denotes 
the full ranges of the 50% completeness level. Stars in the region 
between the two dashed lines are used to determine the edge of 
the TRGB and the peaks of the AGB bump and the RC in the 
following section. 



tion that the / magnitude of the TRGB changes only 
by less than 0.1 ma g for a metallic ity range of —2.2 < 
[Fe/H] < —0.7 dex ()Lee et al.lll993f ). which corresponds 
to the color range defined by the RGB of the metal-poor 
globular cluster NGC 6341 and metal-rich one 47 Tuc, we 
select the stars in this color range (between two dashed 
lines plotted in Figure [9|) to estimate TRGB magnitudes. 
We then calibrate a residual LF, for which contamina- 
tions are removed by comparison with the LF of the 
control field. As deduced from Figure [TOl the final se- 
lected LFs of the GSS's stars show the leading edge of 
the LFs, indicating that the number density of RGB stars 
significantly exceeds the Poisson noise. In fact, the lo- 
cation of the TRGB on the LF is at the point that the 
contamination-subtracted LF with the vertical Poisson 
noise drawn in the middle panel of Fig. [10] is the most 
steeply rising up from the noise. Also, as shown later in 
Figure [HI the global shape of the LFs allows us readily 
to identify the TRGB. 

Table [2] shows the distance to our two stream fields, 
on the assumption that the absolute /-band magnitude 
of the TRGB is M™^^ = -4.1 ± 0.1 for metal-poor 
TRGB stars. This Mj^^^ adopted here corresponds to 
the /-band brightness of the TRGB at [Fe/H] = -1.3 de- 
rived from various empirical and semi-empirical calibra- 



tions CFe rrarese et" ani2000 yfirraro et al.ll20 00':'Lee et al.' 
|l993; Sal aris fc Ca ssisi 199|) bv lSalaris fc"C assisi (2005). 
The TRGB absolute magnitude suffers from statistical 
uncertainty of the order of ~ 0.1 mag in each data. 



Name 


^0,TRGB 


Distance Modulus 


Distance (kpc) 


GSSl 


20.63 ±0.05 


24.73 ±0.11 


883 ± 45 



In estimating the extinction-corrected apparent magni- 
tude of TRGB, /o,TRGB, systematic errors also arise, 
related to zero point uncertainties, calibration through 
relative photometry, aperture corrections, photometric 
errors, smoothing of the LFs and extinction law; the re- 
sulting total error is evaluated by a root-mean-square of 
these errors. 

Our estimated distance to the GSS in the GSSl field 
is 883 ± 45 kpc (see Tab. [2j). It suggests that the 
GSS is located beh ind M31 (which is at 770 kpc by 
iFreedman fc Madara (^199Q) ) and extended along the line 
of sight, whic h is in good agre ement with the distance 
to Field 5 in McConnachie et ajj (|2003i) within errors 
(840 ± 20 kpc), the closest field to GSSl. 



z 10 



n \ ^ ^ ^ ^ \ r 

GSSl 

ITRCB = 20.63 ± 0.03 (m, 




20.5 

Iq (mag) 

Fig. 10. — TRGB detection in the GSS. The top panel presents 
the smoothed LFs as a function of the extinction-corrected /-band 
magnitude zoomed in around the TRGB magnitude. The sohd 
hue shows the original GSSl LF before subtracting the background 
and the dotted line shows the background LF (which is the SE3 
field). The middle panel shows the backgrond-subtracted GSSl 
LF with Poisson error bars. The bottom panel indicates the Sobel 
filter response to the background-subtracted LF. The TRGB mag- 
nitude (the vertical dashed line) including the smoothing error is 
also shown. 



3.3. Metallicity Distributions 

We derive the metallicity distribution (MD) of the 
RGB stars in the GSS, based on the comparison with 
the RGB templa tes defined by Galactic globular clusters 
(|Bellazzini et al.ll2003l) . In this procedure, the assem- 
bly of a sufficient number of the RGB stars with small 
photometric errors is important for deriving an accurate 
MD, otherwise large photometric errors lead to the con- 
tamination of red or metal-rich stars; our Suprime-Cam 
data combined with the careful correction for contamina- 
tions shown in § 12.3.21 are advantageous in this respect. 
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TABLE 3 

Galactic globular cluster data used in the 
present study 



[Fe/H]cG E(B-V) (m - M)o 



NGC 6341^ 


M92 


-2.16 


0.02 


14.74 


NGC 6205^ 


M13 


-1.39 


0.02 


14.38 


NGC 5904^ 


M5 


-1.11 


0.03 


14.31 


NGC 104^ 


47Tuc 


-0.70 


0.04 


13.29 


NGC 6553^ 




-0.16 


0.84 


13.44 


NGC 6528^ 




0.07 


0.62 


14.35 



Note. — The metallicities of these clusters are cah- 
brated in the lCarretta Grattonl (|1997l ) abundance scale. 

^ fSavaane et al.l (l2000l ) 
^ ISagar et all (|T999r) 

^ Locus read off from Fig. 13 of ' Bellazzini et al.l (|2003h 

It is al so remarked that recent studies by iBrown et alJ 
(|2QQ6al lbl) suggest an intermediate age, i.e., younger than 
10 Gyr, for about half the population of the GSS stars, 
in contrast to Galactic globular clusters. However, for 
low-mass RGB stars with age of 6-14 Gyr, the redward 
shift of an RGB sequence with increasing age is minor, 
compared with its sensitive redwar d shift with increas- 
ing metallicity (iGirardi et al.ll2002l). This is a favorable 
behavior of the RGB, lending support to our procedure 
of deriving an MD in comparison with old globular clus- 
ters, but the photometric results should be followed up 
by spectroscopic observations. 

As RGB templates we adop t the six ridge l ines o f 
Galactic globular clusters as in iBellazzini etall (|2QQ3r ). 
The metallicities of these clusters are calibrated in the 
iCarretta fc Grattonl (|1997l ) abundance scale, hereafter 
denoted as [Fe/H]cG- We take four me tal-poor clusters in 
the Galactic halo from Savian e et all (2000): NGC 6341 
([Fe/H]cG = -2.16), NGC 6205 ([Fe/H]cG = -1.39), 
NGC 5904 ([Fe/H]cG = -1.11) and 47 Tuc ([Fe/H]cG = 
—0.70). In addition, we take two metal-rich clusters in 
the Galactic bulge, NGC 6553 ([Fe/H]cG = -0.16) from 
ISagar et all ([l999l) and NGC 6528 ([Fe / H]cG = 0.07, 
ICarretta et al.ll2001f ) fromlOrtolani et alj (|1995f ) (see Ta- 
ble [3|). For secure determinations of MDs using these 
templates, we select the targeted RGB stars having 
-3.8 < Mi < -2.0 and 0.9 < {V - I)o < 4.0 as shown 
in Figure [Til These selection criteria also allow us to 
remove a number of contaminations, such as AGB stars, 
AGB bump stars, young stars, and foreground and/or 
background objects. We perform the interpolation pro- 
cedure to obtain the metallicities of the stars both in the 
stream and control fields, and subtract the derived dis- 
tribution of metallicities in the control field from that in 
the stream field in order to remove the effects of the M31 
field halo stars. 

Figure Ha shows the MDs in the GSSl field of the GSS. 
The vertical error bars denote a nominal uncertainty 
in each metallicity bin, derived from the Poisson errors 
equal to ± Y^A/'(stream + control). It is worth noting that 
these errors are significantly small because of a large 
number of the RGB stars (7051 stars in the GSSl field) 
available from our Suprime-Cam data. It follows that the 
MDs have a broad distribution ranging from [Fe/H] ~ — 3 
to the near solar metallicity and there is a clear high- 
metallicity peak at [Fe/H] peak ^ -0.3. In the GSSl field, 
the average metallicity is [Fe/II]j-Qean — 

-0.55 ±0.45 and 



the median metallicity is [Fe/H]med = —0.45 with a quar- 
tile deviation of 0.23 dex. The average metallicities de- 
rived here are in good agreement with those of kinemati- 
cally selected R GB stars in the GSS (Guhatha kurta etld] 
12006; Kalirai et al. 2006b). 

While the MD is dominated by metal-rich stars, the 
presence of metal-poor stars reaching [Fe/II]< —2 is sug- 
gested from a long metal-poor tail in the MD, based on 
the existence of an extended blue horizontal branch in 
the CMDs derived by HST/ACS (Ferguson et aUEool; 
iBrown ei~ al. 2006a. b). Our found metal-poor tail fea- 
ture is consistent with Bellazzini et al.l ((2003). However, 
as the also discussed, we cannot strongly claim really and 
scattered RGB stars with [Fe/H]> -2. 

3.3.1. Effects of Adopted Assumptions on Metallicities 

The detailed form of the MDs is affected by several as- 
sumptions in the analysis, such as adopted distance mod- 
ulus (m— M)o and reddening correction E{B — V). To see 
these effects in the MDs, we examine the GSSl field as an 
example and the results are shown in Figure [131 While 
Fig. [E] shows the case based on our standard assump- 
tions, the left four panels in Figure [13] show the MDs 
associated with the change of our assumptions. For refer- 
ence, we also plot the MD derived based on the standard 
assumptions as the dotted histogram. Panels (a) and (b) 
are devoted to the effect of changing ±0.11 mag in dis- 
tance modulus, resulting in the change of mean/median 
metallicity at only ±0.07 dex level. Somewhat large 
change of the MDs is seen in panels (c) and (d), where 
we test ±0.05 mag variation in E{B — V). This is in- 
evitable because our method for deriving MDs is based on 
the comparison between the color distributions of RGB 
stars and that of adopted templates, thereby a change 
of E{B — V) affects this comparison in a sensitive man- 
ner, especially in a metal-poor range below [Fe/H] = — 1, 
where the dependence on color is large. The resultant 
change in metallicities is however no more than ±0.18 
dex in response to a ±0.05 mag change in E{B — V). 

Next, we examine the effect of adopting different RGB 
templates, i.e., other than Galactic globular clusters, 
to estimate the metallicity of each star. To do so, we 
re-construct the MDs using the Victoria -Regina theo- 
retical isochrones from IVandenBerg et al.l ([2Q06i) (whose 
isochrones are partly displayed in Fig.[9|), in place of the 
globular cluster templates of the Milky Way. In accor- 
da nce with the interpolation and extrapolation scheme 
of Kalirai et al.l p0 06b). we calculate the metallicity for 
each star in the same segment of the CMD (see also 
their Fig. 6). The resultant MDs with different a- 
enhancement, distance modulus and age are shown in 
the right four panels of Figure [131 Panels (a) and (b) 
are devoted to the effect of changing 4 Gyr in age, re- 
sulting in the change of mean/median metallicity at ap- 
proximately 0.1 dex level. Somewhat large change of the 
MDs is seen in panels (a) and (c) with the two different 
[a/Fe]. The MD with a no-enhancement is very similar 
to the one derived from the Galactic globular cluster tem- 
plates, while the a-enhancement MD has a more metal- 
poor peak. Finally, in panels (a) and (d), we test the ef- 
fect of varying distance modulus. The difference of mean 
metallicity between the two is 0.06 dex. The four MDs 
derived from the model have gradual increase to the peak 
at [Fe/H] - -0.3 without sharp edge at [Fe/H]- -1.1 



A Subaru/Suprime-Cam Survey of M31's Stellar Halo 



15 



GSSl Control 




1 2 3 4 1 2 3 4 



(V - I)„ (mag) (V - l\ (mag) 

Fig. 11. — CMDs of the upper RGB part for two fields. The ridge fines of template Galactic globular clusters are superimposed to each 
plot as solid lines, which are from left to right, NGC 6341 ([Fe/H]= -2.16), NGC 6205 ([Fe/H]= -1.39), NGC 5904 ([Fe/H]= -1.11) and 
NGC 104 ([Fe/H]= -0.70), NGC 6553 ([Fe/H]= -0.16) and NGC 6528 ([Fe/H]= -0.07). The inner square frame encloses the stars whose 
metallicities are determined using the interpolating scheme (jBellazzini et al.ll2003l) . 
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Fig. 12.— Metallicity distribution for the GSSl field. The num- 
ber of stars used to derive it are shown in the upper left corner 
of the panel. The average metallicity with the standard deviation 
and the median metallicity with the quartile deviation of the GSSl 
field are shown in the panel. 

and —0.7 as in MDs constructed based on the Galactic 
globular cluster templates, although all the MDs derived 
from different metallicity indicators have the same global 
shape such as the high metallicity peak at [Fe/H] ~ —0.3, 
the existence of stars with near solar metal and the long 
metal-poor tail. Furthermore, the GSS's mean metallic- 
ity shown in panel (a) of the right of Figure [131 where 
the basic assumptions behind this is almost the same 
as those of Brown et al.l ((2006b), is in good agreement 
with the mean metallicity of ([Fe/H]) = —0.7 derived 
from their best-fit model (based on the distance modu- 
lus of 780 kpc and the Victoria- Regina isochrones). It 
is noteworthy that our mean metallicity estimated from 
only RGB stars is consistent with their result by main- 
sequence fitting to much deeper CMD. 

Based on these experiments, we conclude that the basic 
features of the MDs derived here are robust, including 
the existence of stars with near solar metallicities, high 
metallicity peak at [Fe/H]~ —0.3, and long metal-poor 
tail. The mean and median metallicities remain basically 



unchanged within a typical error of at most 0.1 to 0.2 dex 
in the current method. 

3.4. Stellar Populations Inferred from the AGB Bump 
and Red Clump 

In addition to the broad RGB feature of the GSS, be- 
ing attributable to a composite population of metal-rich 
and metal-poor stars, the CMD shows two noticeable 
populations (see Fig. [9]), the AGB bump and RC, which 
are related to more advanced stellar evolutionary phases. 
To highlight these features, we show, in Figure [HI the 
/-band LFs of the GSSl field, which are corrected for 
extinction. The dashed and solid lines denote LFs be- 
fore and after correction for incompleteness, respectively. 
The vertical dotted line indicates a 50 % completeness 
magnitude of /q = 25.22 for the GSSl field. The main 
features are marked by vertical solid lines: the TRGB, 
AGB bump, and RC. More details on the latter two fea- 
tures are discussed below. 

3.4.1. AGB Bump 

The AGB bump is a stellar evolutionary phase of low- 
mass stars, and is caused by the clustering of stars in a 
CMD, where AGB stars at the beginning of helium shell- 
burning evolution slowly re-ascend along the Hayashi line 
in a thermal nonequilibrium state. However, due to the 
very short lifetime of the AGB phase, the AGB bump 
feature can be d etected only when de aling with a large 
sample of stars (jRejkuba et al.l [20051 ). In fact, iGallartI 
(1998) fir st d e tected an AGB bump from the data of 
^Holland et al.l (|199 6^ in the dense inner part of M31's 
halo (at 7.6 kpc from the center approximately along the 
SE minor axis). As already reported in § 13. H we have 
also found an AGB bump feature in the stream field at 
lo - 23.5 and (V - /)o - 1.2 (Fig. [9]). 

To more quantitatively identify the mean magnitude 
of the AGB bump, we perform a Gaussian- weighted fit 
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[(m-M)o + 0.11] - 7672 stars 
. [Fe/H]„„„=-0.62±0.48 
[Fe/Hi„,j=-0.50±0.26 



[E(B-V)+0.05] - 7521 stars 
. [Fe/H]„„„=-0.73±0.56 
[Fe/H]„,,= -0.58±0.29 \\ 




' ' ' I ' ' ' ' I ' ' ' ' I 
[(m-M)o-O.ll] - 6401 sta^j 
[Fe/H]„^3„=-0.49±0.43 
[Fe/H]„,j=-0.39±0.22 
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[E(B-V)-0.05] - 6506 stat^ 
■ [Fe/H]„„„=-0.42±0.37 
[f"e/H]„,,= -0.33±0.20 
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[a/Fe] = +0.3, dm=24.73, 8Gyr 
. [Fe/H]„„„=-0.68±0.46 . 
[Fe/H]„,j=-0.55±0.24 




[a/Fe] =+0.0, dm=24.73, 8Gyr 
. [Fe/H]„„„=-0.55±0.45 
[Fe/H]„,,= -0.42±0.23 




' ' ' I ' ' ' ' I ' ' ' ' I ' 
[a/Fe] = +0.3, dm = 24.73. 12Gyr 
. [Fe/H]„,^„=-0.78±0.48 
[Fe/H]„,j=-0.66±0.24r 




[a/Fe] =+0.3. dm=24.62. 8Gyr 
- [Fe/H]„„„=-0.62±0.45 
[Fe/H]„,,= -0.49±0.23 
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Fig. 13. — Effect s of varying assumpti ons on the form of the MDs. The dotted histogram: our standard case with (m — M)o = 24.43 and 
E(B — V) from the lSchlegel et all (119981) map, which is derived based on the Galactic globular cluster templates. Left: Panel (a): (m — M)o 
increased by 0.15 mag, (b): (m — M)o decreased by 0.15 mag, (c): E(B — V) increased by 0.05 mag, and (d): E(B — V) decreased by 0.05 
mag. Right: Re-constructed MDs based on the Victoria- Regina theoretical isochrones from VandenBerg et al. (2006 ) with the change of 
[a/Fe], (m — M)o and age. 
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Fig. 14. — /-band luminosity function corrected for extinction 
and incompleteness at magnitude brighter than horizontal branch 
of GSSl field, is plotted with solid line, while the dashed curve 
shows the curve with extinction-correction only. The vertical dot- 
ted lines indicate 50% completeness magnitude of Iq = 25.21 for 
GSSl field. The main features are indicated by vertical solid lines: 
the TRGB at Iq ~ 20.63, the AGE bump at Iq ~ 23.37, and the 
RC at lo - 24.52. 



to the /-band LF of GSS stars after correction for ex- 
tinction and incompleteness. The upper left and right 
panels of Figure [15] shows the zooming CMD and the 
LF with Poisson errors for the GSSl field near the AGB 
bump magnitude, respectively. In this procedure, the 
surrounding stars in the shaded zones (left-hand panel) 
are avoided to reduce unwanted contaminations, since 
the AGB bump feature is rather weak. At the magnitude 
level of the AGB bump, our photometric completeness is 
more than 80 % for /-band imaging. Although, in the 
usual CMD of the left panel, it is somewhat difficult for 
us to detect the AGB bump feature by eye, we can easily 
confirm the existence of the stellar population specified 
by the AGB bump through the LF of the right panel. A 
red solid curve drawn in the LF is the best-fit model with 
the sum of a Gaussian function and a straight line in the 
range of 23.1 ^ /q ^ 23.8. The best-fit function has been 
derived by minimizing reduced chi-squares in the fitting 
procedure. The mean magnitude of the AGB bump is 
then derived from an average value from these fittings: 
for the GSSl field this is /q = 23.37 ± 0.16, correspond- 
ing to M^^^^ = -1.37 with a standard deviation of 0.13 



mag. The estimation of this error is based on those in the 
photometry, smoothing and fitting procedures as well as 
the distance error using the TRGB as shown in § 13. 2[ We 
note that the fitting error in the GSSl field is negligibly 
small (less than 0.01 mag). 

3.4.2. Red Clump 

The RC is a clustering feature of RHB stars being 
met al-rich_and/or_ inter mediate- age. HST/ ACS obser- 
vations (| Brown et al. 2006a. b) of the stream field, which 
is located closer to the galactic center than ours, sug- 
gest the dominance of RHB stars in the HB, i.e., the 
presence of RC. We have also found this feature in our 
stream fields, as already reported in § 13. If at / ~ 24.5 
and (V-/) - 1.0 (Fig. [9]). 

As in the case of the AGB bump, we set tighter limits 
on the mean magnitude of the RC based on a Gaussian- 
weighted fit to the suitably corrected /-band LF. The 
bottom left and right panels of Figure [15] shows the zoom- 
ing CMD and the LF near the RC magnitude, as in the 
upper panel of Figure [15] for the AGB bump. We note 
that particular attention must be paid to the magnitude 
completeness since the RC is just beyond the detection 
limit: the completeness at RC is about 60 % for /-band 
image. We thus reconstruct the completeness-corrected 
LF, using the completeness curves for the stars in the 
unshaded zone of the CMD (the bottom left-hand panel 
of Figure [T5|) . A red solid curve drawn in the LF is the 
best-fit model with the sum of a Gaussian and a straight 
line in the range of 23.8 ^ /q ^ 24.9. Using the curve 
we determine the mean magnitude of the RC: for the 
GSSl field this is Iq = 24.52 ± 0.22, corresponding to 
Mf^^ = —0.21 with the standard deviation of 0.17 mag. 
This error is derived in the same manner as for the AGB 
bump; in this case the photometric uncertainty is larger. 

3.4.3. Age Estimation 

We estimate the mean age of the GSS's stars us- 
ing the relation between age and metallicity for the 
AGB bump and RC features. For thi s purpose, we 
adopt the calibration of this relation by iReikuba et al.l 
(2005) usin^ the theoret ical stellar evolutionary tracks of 
Pietri nferni et al] (|2004l ). Applying to their observations 
of NGC 5128's halo, they showed that the mean age of 
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Fig. 15.— Zoom up display of the CMD of the GSSl field a round 
the A GB bump and the RC features referring to Reikuba et al.l 
(|2Q05h . We do not use the stars in shaded areas to evaluate the 
luminosity functions (see text). Shown on the right panels are the 
/o-band completeness-corrected LFs around the AGB bump and 
the RC, best-fitted with a simple model comprising of a Gaussian 
plus a straight line. The peaks of the AGB bump and the RC 
for the GSSl field are at /q = 23.37 ± 0.16 with a = 0.13 and 
/o = 24.52 ± 0.22 with 0.17, respectively. 



the halo stars in NGC 5128 is estimated as ~ 8^3 5 Gyr. 

Figure [16] shows the absolute /-band magnitudes of the 
AGB bump and RC features predicted by the theoreti- 
cal models, plotted as a function of age for five different 
metallicities. The horizontal dotted lines show the mea- 
sured absolute magnitudes of the AGB bump and RC, 
respectively. If we assume that the alpha elements of 
the GSS's stars are enhanced such that [a/Fe] = +0.3 as 
seen in the local halo stars of the Milky Way, where the 
conversion, [M/H]~[Fe/H]+0.3, holds for most metallic- 
ity ranges, the average metallicity of the stream fields 
is estimated as [M/H] = —0.25, which suggests (as seen 
from asterisks in the Figure) the average age of the GSS 
stars of ~ 8 Gyr. 

To be more quantitative, we compute the average mag- 
nitude of the AGB bump and RC as a function of age, 
as expected fo r the o bserved MDs using equation (9) of 
iRejkuba et al.l (|2005f ) . Figure [T71 shows the correspond- 
ing average / magnitudes of AGB bump and RC features 
plotted as a linear function of age, where the error bars 
are evaluated by employing ±0.1 dex shift in the ob- 
served MDs. The dotted horizontal lines indicate the 
measured mean magnitudes of the AGB bump and RC. 
These magnitude measurements have a negligible statis- 
tical error due to a huge number of available stars, but 
they contain a systematic error associated with redden- 
ing and distance uncertainties, which amounts to about 
±0.22 mag. It is worth noting that for metal-poor and 
old populations with [M/H] < — 1 and age of > 9 Gyr, 
the AGB bump feature is invisible in the theoretical lumi- 
nosity functions, and only a very weak feature is present 
for more metal-rich stars (Re ikuba et al.l [2005 1). so the 
model predictions for the AGB bump with older ages are 




log(age) 

Fig. 16. — Absolute / magnitudes of the AGB bump and RC fea- 
tures predicted by the theoretical models, plotted as a function of 
age for five metallicities (kindly provided by Dr. Marina Rejkuba). 
These metallicities cover most of the metallicity range of the ob- 
served MD in the halo of the GSS. For comparison, the measured 
AGB bump and RC magnitudes are indicated with horizontal dot- 
ted lines. 

not shown. 
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Fig. 17. — Absolute / magnitudes of the AGB bump and RC 
features vs. age predicted by the theoretical models, as obtained by 
convolving the metallicity dependenc y with measured m etallicity 
distributions of the RGB stars (see Reikuba et al.l [20051). Filled 
circles present the relations computed with [a/Fe] = +0.3, while 
open circles indicate those computed with no a-enhancement. The 
measured AGB bump and RC magnitudes are shown as dotted 
horizontal lines. The error bars on the model points are calculated 
from assumed zbO.l dex uncertainty in the measured metallicity 
distributions. 



The mean magnitude of the RC in GSSl field indicates 
that the average age of stars is 7.1 Gyr, whereas the AGB 
bump magnitude in the field suggests younger age. How- 
ever, as seen in Figure [ITJ the dependence of the AGB 
bump magnitude on age is much weaker than that of the 
RC, so that small systematic errors in the measured AGB 
bump magnitude can lead to a large difference in the es- 
timated age. The abundance of a elements somewhat 
affects the AGB bump-based age estimation, in such a 
manner that lower [a/Fe] ratios yield older ages (see open 
circles for [a/Fe]= in Fig. [TTj), whereas the RC-based 
age estimation is rather insensitive to this effect. Much 
larger systematics may be caused by the followings: (1) 
The AGBb feature appears only when [M/H]> —1 and 
age of < 9 Gyr, so the mean age of the stars is dom- 
inated by relatively metal-rich and thus bright popula- 
tions. (2) If the GSS contains relatively young popula- 
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tions associated with secondary star formation activities, 
the mean luminosity of the AGBb reflects the relative 
fraction of each population, especially being biased in fa- 
vor of younger, brighter populations. This effect would 
be more significant for AGBb than RC, because the num- 
ber of the stars in the AGBb feature is much small, so 
the Poisson error in the number of the populations, es- 
pecially for the small number of bright stars among the 
AGBb stars is relatively large, thereby having a larger 
probability of yielding a brighter AGB magnitude and 
younger ages. 

Thus, while taking into account that AGB bump mag- 
nitudes yields younger ages, we adopt the RC-based age 
estimation, giving upper limits on the age of GSS stars. 
Based on the available /-band magnitude of the RC (be- 
cause its F-band magnitude is unavailable in the cur- 
rent study), we arrive at 7 Gyr for the mean age 
of GSS stars, where the error includes 0.22 mag sys- 
tematic uncertainties in the RC magnitude arisen from 
reddening, distance modulus and so on, while the un- 
certainties in stell ar evolut ionary models are not taken 
into account. Reikuba et al. (2005) compared Teramo 
models (Pietr inferni et al. 2004) with Padova models 
(|Girardi et al.|[2000( ) in their similar study for the stellar 
population of NGC 5128 and found an additional 1 Gyr 
uncertainty due to the age difference between the two 
model predictions. 

We note that the mean age of GSS stars derived from 
the RC magnitude is sligh tly you nger compared to that 
derived by B rown et al.i (|2006bl ). They estimated the 
GSS's mean age to be 8.7 Gyr and the mean metallicity 
to be ([Fe/H]) = —0.7 based on the main-sequence fitting 
with the distance modulus of 780 kpc and the Victoria- 
Regina isochrones, whereas our re-computed mean age 
based on the MD constructed by the Victoria- Regina the- 
oretical isochrones (where the MD is shown in the panel 
(a) of the right of Fig. [13]) is 8.1 Gyr. This result is con- 
sistent with the GSS's mean age found bv iBrown et al.l 
(|2006bf ). In addition, the mean age and metallicity esti- 
mations based on different assumptions as discussed in 
§ 13.3.11 are listed in Table [H In particular, all mean 
ages derived from the AGB bump magnitude are much 
younger than those derived from the RC magnitudes. 
The theoretical models provide robust mean ages without 
assumptions such as a-enhancement and age to estimate 
the metallicity, while the ages are somewhat sensitive to 
the assumption of distance as well as in the case of using 
the Galactic globular cluster templates. Therefore, we 
note that our technique applied to calculation of mean 
age based on the RC magnitudes produces several Gyr 
uncertainty in accordance with assumptions such as dis- 
tance, extinction, a-enhancement and stellar evolution- 
ary models. 

The use of the AGB bump magnitude suggests sys- 
tematically younger ages, although this age estimation 
is rather sensitive to the accurate determination of the 
AGB bump magnitude. We note that the presence of 
the AGB bump feature implies the dominance of y oung 
and m etal-rich populations in the GSS; Reikub a et al.l 
(|2005l ) reported that in the models with ages > 10 Gyr 
there exist very few AGB bump stars compared to the 
RC or horizontal-branch stars and that in metal-poor 
([M/H]< -1 dex) and old (age > 9 Gyr) models the 



AGB bump is not present at all. Thus, the presence of 
the AGB bump itself yields the younger age estimate. 
It is also noted that the current implication for young 
and metal-rich populati ons in the GSS m ay be associ- 
ated with the result of Brown et al." (2006b') (as depicted 
in their Figure 17b), suggesting that star formation ac- 
tivities in the GSS had continued until some Gyrs ago; it 
is inferred that the progenitor of the GSS may have had 
three main star formation activities: (1) 30% of the stars 
were provided in the first star formation (e.g., by the col- 
lapse of the initial gas) before 8 Gyrs ago, (2) 40% of the 
stars were provided in the second star formation (e.g., su- 
pernova feedback or galactic interactions) about 8 Gyrs 
ago, and (3) remaining 30% of the stars were provided 
in the third star formation (e.g., the tidal interaction in 
the first pericentric encounter with M31) after the second 
star formation. 

3.5. Mass of the Progenitor Galaxy 

The total mass of the stellar system expected for the 
progenitor galaxy of the GSS can be estimated from 
the mean metallicity of stars in the GSS ([Fe/H] mean = 
—0.7), combined with the metallicity- lumin osity relation 
for Local Group dwarfs (jCote et al.l l2000l ). The rela- 
tion shows, on average, more metal-rich populations in 
more luminous galaxies, and we obtain Mgtars ~ 10^ 
M0 for the progenitor galaxy if the mass-to-light ratio 
is M/L ~ 2. In contrast, assuming the stellar mass- 
met allicityre]^;ti^n_obseryed^ the Local Group dwarf 
galaxies of lDekel fc Wool ([2003h . we find M - 4 x lO^M© 
as the prog:enitor ni a ss, w hich is consistent with the 
work of iMori fc RichI (|2008l ) who restricted the progen- 
itor mass to be M < 5 X lO^M©, taking into account 
the effect of disk heating by dynamical friction. Such a 
dwarf galaxy may have been accreted to the M31 halo at 
recent epochs, possibly within some billion years, leaving 
yet a conspicuous giant stream like the Sgr stream in the 
Milky Way. 

4. PANORAMIC VIEWS OF THE ANDROMEDA STELLAR 

HALO 

In this section, we report on the results for the entire 
halo regions of M31 observed by the current survey, i.e., 
including those along its minor axis both in the north 
and south parts of the galaxy and a south-west major- 
axis region. 

4.1. The Color- Magnitude Diagram 

Figure [18] shows a co-added CMD from all of the 
observed fields (including the GSS field as shown in 
the previous section), after removing extended sources 
such as background galaxies an d cos mic rays based on 
DAOPHOT parameters (see § [2X21) . The thin solid 
lines in the CMD show theoretical RGB tracks from 
IVandenBerg et all (2006) for an age of 12 Gyr, [a/Fe] = 
+0.3, and metallicities {left to right) of [Fe/II]= —2.31, 
-1.71, -1.14, -0.71, -0.30 and 0.00. The dotted line 
denotes about 90 % completeness limit as determined 
by artificial star experiments, where the limit is esti- 
mated by making use of the NW8 field with worse see- 
ing and shorter exposures. Although redundant stars 
and galaxies are not completely removed from the CMD, 
we can clearly identify a broadly distributed RGB fea- 
ture attributable to overdense regions with a number 
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TABLE 4 

Mean age and metallicity based on various assumptions 



Various Assumption ([Fe/H]> {age>Rc (age>AGBb 



Based on the Galactic globular cluster templates 



Standard 


-0.55 


7.1 


0.7 


Standard with no a-enhancement 


-0.55 


8.2 


2.3 


(m — 


M)o +0.11 (24.84) 


-0.62 


4.8 


1.0 


(m — 


M)o - 0.11 (24.62) 


-0.49 


9.4 


2.8 


E(B - 


-y) + 0.05 (0.10) 


-0.73 


5.6 


0.4 


E{B - 


-V)- 0.05 (0.00) 


-0.42 


8.9 


1.4 






Based on the Victoria-Reg 


;ina theoretical isochrones 






a/Fe 


= +0.3, dm= 24.73, 8 Gyr 


-0.68 


8.0 


1.9 




a/Fe 


= +0.3, dm= 24.73, 12 Gyr 


-0.79 


8.1 


2.3 




a/Fe 


= +0.0, dm= 24.73, 8 Gyr 


-0.55 


8.2 


2.5 




a/Fe 


= +0.3, dm= 24.62, 8 Gyr 


-0.62 


9.4 


4.9 



of metal-rich stars such as the GSS. At fainter than 
Iq ~ 23 (mag), the RGB feature appears to be much 
broader because of the presence of background galaxies 
which passed through our selection of separate stellar 
objects. Although it is contaminated by these objects, 
a population extending toward fainter and bluer region 
than the most metal-poor theoretical RGB isochrone at 
Jo ^ 24 (mag) corresponds to the horizontal branch stars 
in M31's halo. 

In the right panel of Fig. [18] which is the same CMD 
as the left one but with the higher maximum level of 
representation in the stellar densities, there exists a se- 
quence of Galactic disk dwarfs at (V — /)o ^ 2 having 
a broad RGB; the vertical sequence is the result of low- 
mass stars accumulating in a narrow color range, yet be- 
ing seen over a large range in distance along the line of 
sight. In addition, halo stars in the Milky Way can be 
seen as vertically-distributed stars on the blue side of 
this diagram, (V — I)o < 0.8 and Iq < 23. Usually this 
appears as a smooth vertical structure in a CMD, which 
corresponds to the stars at or close to the main-sequence 
turnoff at increasin g distance through the Galactic halo 
(jMartin et al.ll2QQ7l : Tlbata et al.ll2QQ7[ ). 

4.1.1. Calibration of the Control Field 

Our three control fields, to be utilized for removing 
contaminations such as Galactic dwarf stars which can- 
not be separated morphologically, are located at the same 
galactic longitude / = 103° but at different galactic lat- 
itudes with b = -16?5, -20?7, and -25?1, i.e., far from 
the M31 direction (/ = 121.2°, 6 = -21.6°). The CMDs 
of the three control fields are shown in Figure [191 from 
left to right, galactic latitude is higher. For reference, 
the solid/dotted white lines as shown in Fig {18] are also 
plotted in these CMDs. A prominent feature is that the 
number of Galactic disk dwarf stars is gradually decreas- 
ing with increasing galactic latitudes. In contrast, there 
is little change in the distribution of remaining faint back- 
ground galaxies between these three CMDs. Finally, it 
is found that the number variation of the Galactic halo 
stars in these CMDs is also small. 

These properties of the control fields indicate that the 
contaminants in our target fields of M31's halo are dom- 
inated by the faint background galaxies and the Galac- 
tic disk dwarf stars, while the effect of the Galactic halo 



stars can easily be corrected because such stars distribute 
outside the M31 RGB sequence in the CMD. Regarding 
the unresolved background galaxies, one can assume that 
such galaxies are uniformly distributed in a Suprime- 
Cam field-of-view, whereby these background contami- 
nants per field are simply regarded as having a uniform 
density, regardless of the galactic latitude. Thus, the 
remaining task is to understand the effects of the fore- 
ground Galactic disk dwarf stars on our target fields. 

In order to correct the effect of these dwarf stars in 
each target field of M31's halo, we estimate the fraction 
of such stars in each target field based on the knowledge 
of the number of the contamination stars in the CMDs 
of the control fields, by combining the observed CMD s 
and those predicted by the models of lRobin et al.l (|2QQ3l ). 
Firstly, we use two CMDs of the short exposed NWl 
field and one of the three control fields with almost the 
same galactic latitude of b ~ —21°. We have obtained 
0.78 as the scale factor to adjust the number difference 
arising from the galactic longitudinal difference between 
the M31 with / ^ 121° and our control fields with / ^ 
103°, based on the comparison of the star counts of clean 
Galactic dwarfs in a box defined as 17.0 < Iq < 18.0 
and 1.5 < (^ — /)o < 2.2. Note that these selection 
criteria enable us to avoid the effects of bright AGB and 
RGB stars of M31. This scale factor derived from our 
observation i s cons istent with the model prediction from 
iRobin et al.l (|2003i ) including Poisson noise. In addition, 
the scale factor calculated from the star count of the 
disk stars in a box defined as 20.5 < Iq < 22.5 and 
1.2 < {V - I)o < 3.0 of the corresponding model CMDs 
without the M31 population is also in good agreement 
with the scale factor from the comparison of the bright 
disk stars. Secondly, we investigate the variation of the 
number density of the disk dwarf stars along galactic 
latitude. Figure [20] shows the relative number variation 
normalized by the number at the control field with b = 
—20? 7 as a function of the galactic latitude, where the 
disk stars reside in a box defined as 20.5 < Iq < 22.5 
and 1.2 < {V — I)o < 3.0. The three big solid circles 
denote our observed results, while the small solid circles 
show the model prediction by Robin et al. (2003) with 
the same galactic coordinates as our M31 target fields 
and the solid triangles are the scaled plots by the above- 
mentioned scale factor of 0.78. It is evident that the 
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Fig. 18. — A composite CMD in the Iq vs. {V — I)o plane for stellar sources detected in all the observed fields (including the GSS field). 
The dotted lines denote the full range of the 90% completeness levels for the mosaic images of the Suprime-Cam in the NWS field. The 
superposed solid lines in the CMD show theoretical RGB tracks from VandenBerg et al. (2006) for an age of 12 Gyr, [Q;/Fe]= +0.3, and 
metallicities (left to right) of [Fe/H]= —2.31, —1.71, —1.14, —0.71, —0.30 and 0.00. Extended background galaxies have been removed via 
image classification. The right panel is the same CMD as the left panel, but with the higher maximum level of representation in the stellar 
densities. 
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Fig. 19. — CMDs of the three control fields in the Iq vs. (V — I)o plane. For reference, the white solid/dotted lines as shown in Fig. 1181 
are also overlaid, left: (l,b) = (103?0, -16?5), middle: (l,b) = (103?4, -20?7), right: {l,b) = (103?4, -25?1). 



spatial variations of both the observational result and 
the model prediction are reasonably consistent with each 
other. Based on the results of these exp eriments, we 
thus a dopt the theoretical galactic model of iRobin et alJ 
(|2QQ3[ ) to estimate the spatial variation of the number of 
the disk dwarf stars. This spatial variation of the number 
density is well represented by an exponential profile at 
the lower galactic latitude than that of the NWl field and 
is represented by a linear profile at the higher galactic 
latitude, as judged by a (dotted) fitting line in the figure. 
For this correction of the contaminations, we adopt the 
blank field with the lowest galactic latitude of 6 = —16? 5 
as the control field for the north part of the M31 halo, 
whereas the blank field with the galactic latitude of b = 
—20? 7 is adopted as the control field for the south part 



of the M31 halo. 

4.2. Stellar Population Maps 

In this subsection, we investigate how the RGB pop- 
ulations of the M31 halo, as identified in Fig. [181 ^ire 
distributed spatially out to about 100 kpc from the M31 
center. To do so, we extract these stellar populations of 
M31 alone from our imaging data. However, it is noted 
that the outer halo of M31 is known to have a very diffuse 
stellar densit y as deduced from pre vious studies (e.g., 
[Gilbert et aH 120061 : llrwin et al.ll2Q05f ). so a clean signal 
of these M31 populations is buried under heavy contam- 
ination from foreground disk dwarf stars and the back- 
ground galaxies. In particular, the amount of the disk 
dwarfs in the northern part of the M31 halo are about 
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Fig. 20. — The variation of the number density of the Galactic 
disk stars as normalized by that in the control field with b = —20? 7 
as a function of the galactic latitude. The disk stars are taken 
from a box defined as 20.5 < Iq < 22.5 and 1.2 < {V - I)o < 3.0 
in the CMD. The three big solid circles denote our observational 
results, while the small solid circles show the model prediction by 
iRobin et al. (2003) with the same galactic coordinates as our M31 
target fields and the solid triangles are the scaled plots. The solid 
triangles are fitted by exponential (at the lower latitude than that 
of the NWl field) and linear profiles (at the higher latitude). 

twice as much as in the southern region, as suggested in 
Fig. [20l In order to reveal the signal from these heav- 
ily contaminated data, we adopt th e so-called matche d 
filter method as has been utilized in llbat a et alJ (|2QQ7l ). 
This method is an optimal search strategy if one has a 
precise idea of the properties of both the signal and the 
contamination. 

We first divide the /q vs. {V — I)o plane into uniformly- 
spaced rectangular grids. Then, based on the simple 
weighting of each CMD in a box defined by these grids, 
in terms of the ratio of signal to contamination as ex- 
pected for the concerned CMD, we can optimally boost 
the signal and suppress the contamination. Since we 
have the objective CMDs (i.e., containing both the sig- 
nal and contamination) as well as the CMDs containing 
only the contamination, it is possible to obtain the so- 
called weight matrix, which is given by subtracting the 
CMD for the contamination from the objective CMD (us- 
ing the correction method for the contamination as dis- 
cussed in the previous subsection), and finally we sum up 
these corrected CMDs in each field of M31's halo. The 
weight matrix obtained by this procedure is shown in 
Figu re [22l Th is is reasonably similar to the one derived 
by Ib ata et a l. (2007). Here, to construct the weight ma- 
trix, we have selected the stars filed with 90% complete- 
ness limit, the fainter stars than an upper dashed line of 
Fig. HH] and the stars in the possible metallicity range of 
-2.31 < [Fe/H] < H-0.00, as an example. Note that, as 
expected, the greatest weight arises at faint magnitudes 
in the color range of 0.9 < {V — I)o < 1.5, so of course, 
stars with this photometric property will contribute most 
strongly to the following matched filter maps. 

Based on this method, we have analyzed all of the 
CMDs and obtained the matched filter maps. In Fig- 
ure [211 we present the final matched filter maps over the 
entire survey region for four different ranges in metal- 



licity. The maps are smoothed with a Gaussian kernel 
over 0.054°. The limiting magnitudes are chosen to be 
Jo = 24.0 and Vq = 24.8; these magnitudes keep the de- 
tection efficiency over 90% completeness in the entire re- 
gion. The images with the resolution of 0.018° x 0.018° 
are displayed in a logarithmic scale for the left panels, 
whereas the right panels are shown in a linear scale. The 
inner ellipse overlaid in these maps represents a disk of 
inclination 77° and radius 2° (27 kpc), the approximate 
end of the regular Hi disk. The outer ellipse shows a 
4° (55 kpc) radius ellipse which is flattened to have the 
axis ratio of c/a = 0.6 [i.e ., standard guide ellipse firstly 
drawn bv [Ferguson et al.] (|2002 )]. In addition, the cyan 
concentric dashed circles show projected radii of 80 kpc. 
We note that the logarithmic-scaled maps are useful for 
elucidating the global spatial variation of the M31 pop- 
ulations in the inner halo and the large-scale smooth 
spread of the diffuse outer halo, while the linear-scale 
maps are useful for the presentation of some faint sub- 
structures in the outer halo, with signal boosted by the 
adopted CMD filter. We summarize the characteristic 
features in these maps as follows. 

• South Field: 

The southern quadra nt of the M3 1 halo has already 
been mapped out bv llbata et al.l (|2007l ). We have 
confirmed their finding of some diffuse substruc- 
tures in this part of the halo as given below. 

1. The GSS. This is the most conspicuous over- 
de nse region in the M 31 halo firstly detected 
by Ibata et alj (|2001l ) . As discussed in the pre- 
vious section, the GSS has high mean metal- 
licity of [Fe/H] = — 0.7±0.5 and slightly young 
population with intermediate age of about 8 
Gyr. It is interesting to note that in Fig. [2Th 
for the most metal-poor range, the GSS has 
almost disappeared. 

2. Stream (7 at (^ ~ 3.5,7^ ~ —2.5), which 
was discove red by the MegaCam survey of 
llbata et all (2007). This stream is visible in 
the more metal-rich range except for Fig.[2Th. 

3. Stream D at (^ y 2.5, ?7 ^ —1.8) , which 
was also found by llbata et al.l (|2007[ ) for the 
first time. This strong, stream-like structure 
is detected as a clumpy substructure in the 
intermediate metallicity range (in particular, 
Fig.EIb). 

4. Major-Axis Diffuse Structure at (^ ^ 
—2.5,7^ ~ —4). A subtle overdense structure 
may be a part of the diffuse structure extend- 
ing on the major axis betw een a projected dis - 
tance of 50 and 100 kpc, as llbata et al.l (|2007l ) 
have mentioned in their paper. One can more 
clearly identify the faint substructure in the 
map over the full metallicity range of Fig.[2Tli. 

• North Field: 

Previous works have investigated only the limited 
northern parts of the M31 halo near the center. 
Here, we first explore the regions far from the cen- 
ter, i.e., the outer halo. 
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Fig. 21. — Matched filter maps to a limiting magnitude of Iq = 24.0 and Vq = 24.8. Log-scaled images are shown on the left, while 
linear-scale versions are presented on the right. The resolution of this map is 0.018° x 0.018° pixels, smoothed with a Gaussian kernel over 
3 pixels. 



1. Western Shelf (W-shelf), which is a faint 
stream- hke feature on the western side of M31, 
located at ^ —1.5,?^ ~ 1). This struc- 
ture has already been detected in the map of 
iFerguson et alJ (2002) , having similar color to 
the GSS. There is a sharp cutoff of its stellar 
distribution in the most metal-rich range of 
Fig. [2Bi. 

2. A previously unknown stream is seen at ~ 
—3.5,7^ ~ 3). This structure is most clearly 
identified in the map of the intermediate 
metallicity range, i.e., in Fig. [2Tb . We refer 
to this stream as "Stream E" in the discus- 
sion below. 

3. A further faint low- surface-brightness over- 
dense structure is detected toward ~ 
—5.3, T] ~ 4.3). The increased sensitivity with 
the full metallicity range as shown in Fig. [2Tl i 
has revealed this structure, which we refer to 



as "Stream F" . 

4. The extension of the underlying halo in the 
metal-poor range of the North- West minor 
axis field is comparable to that found in the 
South-East minor axis. This is not the case in 
the metal-rich range as shown in Fig. [2Tk . 

The maps displayed in Figure [21] show the distribu- 
tion of the matched filter statistics, whereby the resulting 
counts are therefore rather difficult to interpret directly. 
This is principally because the matched filter method re- 
lies on a model of the stellar population that one desires 
to detect and the statistics we measure depend on a com- 
bined CMD of our survey fields, which however do not 
cover the entire region of the M31 halo. In the next sub- 
section, we discuss in more detail the populations that 
are highlighted in Figure [211 by determining their funda- 
mental properties such as the metallicity and age based 
on the methods we have already described in the previous 
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Fig. 22. — Matched filter weight map on a log scale, trimmed 
to the color-magnitude region encompassing stars of metallicity 
-2.31 < [Fe/H] < +0.00. The resolution of this map is 0.2 (mag) 
X 0.25 (mag). 

subsection. 

5. SPATIAL SUBSTRUCTURES 
5.1. Western Shelf 

This W-shelf substructure, which has low surface 
brightness compared with even the core of the GSS, was 
already mapped out as the ove rdens e re^on of the RGB 
stars by 'Fergus on etHI ((2002)) and llbata e t alJ ("20071). 
Some A/'-body simulations suggest that several other sub- 
structures in the M31 halo including the W-shelf can be 
explained as the forward continuatio n of the GSS (e.g., 
iFardal et al.ll2007l : iMori fc Richll2008f ). However, despite 
the results obtained by such simulations, the detailed ob- 
servational studies of t he faint W-shelf substr ucture have 
been little made (e.g., Richard son et al.]l2009l ). especially 
to set useful limits on its origin. In this subsection, we 
discuss the fundamental properties of this shelf by com- 
paring them with those of the GSS based on our results 
(§[3]). In the analysis we describe below, it is possible to 
remove systematic errors in this comparison between the 
W-shelf structure and GSS, since the basic data of both 
fields, such as distance, metallicity and age, are extracted 
based on the same technique. 

Figure [23] shows a background/foreground subtracted 
CMD for the W-shelf, combining the catalogs of the 
17 kpc to 24 kpc fields along the North- West minor 
axis. This range is chosen to avoid the contamination of 
NGC 205 and to be outside of the W-shelf edge of - 1.8° 
from the center. This edg e feature w as also reported in 
Figure 1 and Table 1 of F ardal et al.l ([2007|). The super- 
posed solid lines in the CMD show the same theoretical 
RGB tracks as Fig. [181 This CMD has similar morphol- 
ogy to that of the GSS, as represented by a broad/metal- 
rich RGB and a slightly young population attributable 
to a RC extending toward brighter magnitude (see also 
Fig. [8]). However, it is noteworthy that the population 
of the W-shelf with [Fe/H] ~ —0.3 is somewhat sparser 
than the GSS. Furthermore, because of the high quality 
data of the W-shelf compared with that of the GSS, an 
old/metal-poor blue horizontal branch (HB) extending 



toward (V — /)o ~ and Iq ~ 25.2 is clearly identified, as 
well as the brightest part of the RGB bump, which indi- 
cates the presence of metal-rich populations. The slightly 
diffuse feature of the blue HB is simply due to the large 
photometric error of ±0.2 in color. We cannot detect a 
prominent AGB bump as was seen in Fig.[8l probably be- 
cause the W-shelf field has lower surface brightness than 
the GSS field. Note that our CMD of the W-shelf con- 
tains yet an extra contamination of the underlying halo 
component, although the effect appears to be not signif- 
icant if we consider the approximately 1 mag arcsec"^ 
fainter correction for the underlying halo than that for 
the W-shelf (see below Fig. [37|) . 

The distance to the W-shelf has been estimated by the 
sobel edge detection algorithm as was already described 
in § 13. 2[ We have found the distance modulus of 24.51 
mag corresponding to 798 (±40) kpc, on the assump- 
tion that the absolute /-band magnitude of the TRGB 
is M™^^ = -4.1 ± 0.1 for metal-poor TRGB stars. It 
suggests that the W-shelf region probed by our Suprime- 
Cam pointing lies about 85 (±10) kpc relatively in front 
of the GSS probed by our observation. The error in this 
estimate depends only on a small photometric calibration 
error, without considering any systematic errors such as 
the determination of a TRGB absolute magnitude. 

Based on our estimated distance to the W-shelf, 
we now proceed to construct its MD, assuming stel- 
lar population with [a/Fe]= ±0.3 and age of 10 Gyr 
based on the Victoria- Regina theoretical isochrones from 
Vanden Berg et al. (2006). Figure [24l shows the MD and 
its cumulative distributions in the W-shelf field plotted 
as the dashed histogram. The foreground/background 
contaminations have been subtracted based on the scaled 
control field as described earlier. In contrast, solid and 
dotted histograms denote the MDs of the GSS. The dot- 
ted histogram is the same MD as the top-left histogram 
in the right panel of Fig. [131 ^ind the solid histogram 
is the scaled MD of the GSS; the contamination for the 
former dotted MD is estimated based on the SE3 field 
which is located near our observed GSS along the minor 
axis (see Fig. [T]), while the contamination for the latter 
solid MD is derived from a control field far from the M31 
center, that is, free from the population of the underlying 
halo in M31. This consistent correction procedure for the 
contamination makes it possible to deduce the nature of 
the W-shelf itself, based on the comparison between the 
dashed and solid histograms. 

As is evident, these two MDs are basically the same in 
the metal-poor range of [Fe/H] ^ — 1, characterized by a 
prominent metal-poor tail. A blue HB feature identified 
in the CMD of Fig. [23l also suggests that this metal- 
poor tail is attributable not to the contamination of the 
metal-rich AGB stars but to the existence of the remark- 
ably metal-poor and old population. This suggests that a 
progenitor (s) of such overdense substructures in the inner 
part of the halo had undergone primordial star formation 
before it (they) intera cted dynamically w ith M31 2 Gyr 
ago, as simulated by Fardal et al.l (|2007l ). In contrast, 
the properties of these two MDs in the metal-rich range 
appear slightly different, although mean/median/peak 
metallicities of the W-shelf are in agreement with those 
of the GSS, [Fe/H]mean = -0.75, [Fe/H]med = -0.63 and 
[Fe/H] peak = —0.5. Th e two-sided Kolmogorov-Smirnov 
(K-S) test (jPress et al.|[l992.) actually rejects the null hy- 
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pothesis that the two data sets are drawn from the same 
distribution, with a probability less than P = 0.03, which 
is derived from the maximum value of a vertical devia- 
tion between the two cumulative distributions (as in a 
bottom panel of Fig. [24|) . This indicates that these two 
MDs are different in a statistically significant manner. 

We proceed to investigate the mean age of the W-shelf 
population based on the constructed MD and the RC 
magnitude, through the same manner as the analysis of 
the GSS in the previous section. By least-square fitting 
to the LF model with a Gaussian function and a straight 
line to the measured luminosity function, we estimate 
I^^ = 24.31 ± 0.16 and Mf ^ = -0.20, whereby we find 
8.2^3^ Gyr for the mean age of the W-shelf stars. 

Based on these analyses, we can conclude that the fun- 
damental properties of the W-shelf population such as 
the mean metallicity and the mean age are excellently 
consistent with those of the GSS probed by our Suprime- 
Cam pointing. Also, there exists a small but statisti- 
cally meaningful difference in the fraction of metal-rich 
stars between the W-shelf and the GSS. We emphasize 
that this difference in the metal-rich population is ro- 
bust because our data are taken with higher detection 
efficiency, much smaller photometric errors and larger 
statistics than previous studies. To understand the ori- 
gin of this small difference in the metal-rich population, 
it is useful t o consider the simulation results, such as 
iFardal et all (^2007) , which show that the W-shelf struc- 
ture was formed by the third pericentric passage of the 
same progenitor as that produ c ing th e GSS. In particu- 
lar, the model bv IFardal et al.l (|2008l ) suggests the exis- 
tence of a population gradient in the progenitor galaxy, 
i.e., spatial gradient in metallicities, as found mainly in 
disk galaxies. The observational evidence for the po pu- 
lation gradient was discovered by llbata et aH (|2007l ) , in 
their comparison between the core and envelope popu- 
lation of the GSS. Our finding of the small difference in 
the fraction of metal-rich stars between the W-shelf and 
the GSS may be caused by such a population gradient 
within the progenitor galaxy. Of course, there remains 
yet another possibility that the origins of the GSS and 
the W-shelf are different. 

5.2. South-East Streams C and D 

Figure [251 shows close-up maps of the south-east minor- 
axis region in the proximity of our GSS field at ~ 
0.5,77 ~ —2.2), divided [Fe/H] into two non-overlapping 
ranges (+0.00 to -0.71 and -0.71 to -1.71). This 
division for metallicities makes it clear that two over- 
dense regions are indeed present, contrary to the ex- 
pected smoothly declining stellar density of the halo. 
These over dense regions, which were first discovered by 
Ifbata et all (|2007[ ). are called Stream 'D' (at ^ 2.5, rj 
— 1.8) and 'C' (at ^ ~ 3.5, r] ~ —2.5) in order of increas- 
ing declination. Actually, they are aligned in a direction 
perpendicular to the minor axis and more extending like 
a stream feature on both sides of our observed regions; 
we see only a part of the two streams. 

Unfortunately these substructures show about 1/10 
lower surface brightness than the GSS and the W-shelf, 
and the detection limits for these data are somewhat 
worse than for the control field data. Thus, the method 
for subtracting the background contaminations in the 
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Fig. 23. — A CMD in the Iq vs. {V -I)o plane for stellar sources 
detected in the W-shelf field. The dotted lines denote the full 
range of the 50% completeness level. The superposed solid lines 
in the CMD show theoretical RGB tracks from VandenBerg et al.j 
(:2006) for an age of 12 Gyr, [Q;/Fe]= +0.3, and metallicities (left 
to right) of [Fe/H]= -2.31, -1.71, -1.14, -0.71, -0.30 and 0.00. 
The foreground/background contaminations have been subtracted 
based on the scaled control field of CompL103Bml6_5. 
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Fig. 24. — Upper: Metallicity distribution functions (with er- 
ror bars denoting Poisson uncertainties) for the W-shelf sample 
and the GSS sample, as interpolated from the chosen Victoria- 
Regina isochrones. The same photometric limits as Fig. 1121 have 
been imposed. The control field, normalized with the gradient 
shown in Fig. [20] have been used to subtract the expected fore- 
ground/background counts in each of the metallicity bins. The 
dotted histogram is the same MD as the top-left histogram in the 
right-hand panel of Fig. 1131 Lower: The cumulative MDs for the 
two substructures. The solid and dashed distributions appear very 
similar but the K-S statistics suggests that they are different. 
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CMD would need caution, when we simply compare with 
the control field having much better detection efficiency. 
Therefore, in Figure [26] we display the contamination- 
uncorrected CMDs of the stream-like structures within 
the range of projected R ~ 50-60 (35-45) kpc for 
Stream C (Stream D). These CMDs are smoothed with a 
Gaussian kernel over 3 pixels {d{V — I) = 0.03, dl = 0.03) 
in order to diminish the effect of spiky noises attributable 
to low statistics. Here, as a further step, we have tried 
to roughly remove background galaxies in combination 
with a criterion based on the FWHM parameter derived 
from SExtractor (jBer tin & A rnoutsI 11996). The crite- 
rion was chosen along the stell ar sequence in the space 
of magnitude vs. FWHM fe.g.. iKashikawa et all (2004). 
Since the detection method of sources in SExtractor suf- 
fers from blending objects, the detection limit gets into 
somewhat shallower in response to changes in worse see- 
ing size. Nonetheless, both of the CMDs of Stream C and 
D as in Fig. [26] imposed by the stringent criterion present 
a clear RGB sequence along the theoretical isochrones at 
lo < 23.5 and (V - /)o - 1. 

A prominent feature of the CMDs is that Stream C 
is predominantly metal-rich, being comprised of the 
stars in the range -1.14 < [Fe/H] < -0.30, while 
in Stream D metal-poor stars with the metallicity of 
-1.71 < [Fe/H] < -0.71 are weh populated. The 
corresponding MDs of the two streams given in Fig- 
ure [27] statistically confirm these features, where we have 
made the subtraction of contaminations as was done 
for the GSS and W-shelf fields. The mean metallicity 
with the standard deviation of Stream C (Stream D) is 
[Fe/H] = -0.89 ± 0.55 (-1.03 ± 0.61) based on 1684 
(1288) stars, assuming an age of 1 2 Gyr, [a/Fe]= -hO.3 
and distance modulus of 24.43 from Freedman & Madore 
|l990). The characteristics of this population differ- 
ence are also suggested by th e previous study based on 
the CFHT/MegaCam survey (jibata et al.ll2007[ ): we have 
confirmed the existence of these faint outer streams based 
on our Subaru/Suprime-Cam data with higher S/N. 

Now we further determine the mean age of the domi- 
nant populations based on the MD and RC magnitude. 
We present the luminosity functions of the Stream C and 
D in the expected magnitude range of the RC when as- 
suming the aforementioned distance in Figure [28] The 
RC magnitudes of Stream C and D are /q = 24.37 and 
24.47, respectively. By determining a Gaussian peak as 
RC magnitude and convolving the MDs of Figl27] with 
the magnitude, we estimate the mean ages of the dom- 
inant populations of these streams as ~ 9.3l?:^ Gyr for 

Stream C and ~ 9.6t?;{ Gyr for Stream D. The errors are 
based on the photometric errors, which strongly depend 
on the relation between RC magnitude and metallicity 
distribution at younger part as shown in Fig. [16] 

5.3. Major-Axis Structure 

The faint diffuse structure extended along the south- 
west major axis between a projected distance of 50 and 
100 kpc was mentioned in Ibata et al. (2007). Although 
our Suprime-Cam survey regions are slightly different 
from those covered by CFHT/MegaCam and both of 
the two surveys do not encompass the entire region of 
the diffuse major-axis structure, the CMD of these areas 
shows a clearly well-populated metal-poor RGB sequence 



with a dominant population in the metallicity range of 
-2.31 < [Fe/H] < -1.14 (Figure [29|). The correspond- 
ing MD and cumulative MD using about 420 stars and 
the Victoria-Regina theoretical isochrones for an age of 
12 Gyr and [a/Fe]= +0.3 are displayed in Figure [30] 
They are constructed by imposing further color criterion 
with (V — /)o < 1.8 in order to avoid the heavy fore- 
ground disk dwarf contamination, considering the exis- 
tence of the dominantly populated metal-poor stars as 
shown in the CMD. The mean metallicity of the popula- 
tion is [Fe/H] ~ —1.2. Our confirmed stellar population 
in the diffuse structure on the south-west m ajor- axis is 
in nice agreement with that of llbata et al.l ( ^2007 ), and 
there is no population gradient of the diffuse structure in 
b oth of the su rvey fields. 

llbata et all |2007) concluded that the two stellar pop- 
ulations of the diffuse structure and the GSS are dif- 
ferent and likely unrelated despite the visual impression 
that the density map of their Figure 21 gives indications 
that the major- axis population merges with the GSS. 
However, in the LF shown in Figure [HJ we detect a 
subtle peak similar to the RC over Poisson noise, sug- 
gesting that there is an intermediate-age population in 
the diffuse halo structure. If the peak morphology with 
Iq = 24.34 is actually attributable to the RC, we can es- 
timate the mean age of the population as ~ 8.6l5"5 Gyr. 
This estimation provides a comparable age to th e GSS . 
Based on their numerical simulation, Fardal et al.l (j2008l ) 
actually reproduced the decline in the mean metallicity 
from the central core of the GSS to its cocoon to the 
southwest, using an accretion model of a disk-like pro- 
genitor having a strong radial metallicity gradient. Thus, 
this may suggest that the metal-poor /young population 
of the diffuse halo structure originated from the debris 
of the GSS's progenitor. So we cannot readily rule out 
the association of the diffuse structure with the GSS. 
A much deeper and wider survey reaching down to the 
main-sequence turn-off by a planned thirty meter tele- 
scope is necessary to more accurately determine the age 
and the distance to the diffuse structure. Additionally 
the spectroscopic survey to measure kinematics will be 
needed to understand the origin of these halo structures 
in detail. 

5.4. New North- West Streams E and F 

Here, we report on two new substructures seen in the 
north-west portion of the outer halo of M31. In con- 
trast to the south quadrant of M31's halo which has been 
studied by a few other groups, the north quadrant has 
remained unexplored until now, in part because its low 
Galactic latitude {h ~ —10°) results in a high degree of 
contamination by foreground Milky Way dwarf stars. 

The two overdense structures on the north-west mi- 
nor axis at ~ 60 and 100 kpc can be deduced from 
the visual inspection of the matched filter maps in Fig- 
ure [32] This figure shows close-up maps of the north- 
west minor-axis region, divided by [Fe/H] into two non- 
overlapping ranges (+0.00 to -0.71 and -0.71 to -1.71). 
Here, we refer to these as Stream E and F, corresponding 
to the inner 60 kpc and outer 100 kpc substructures in 
the discussion below. Their dominant RGB population 
has a metallicity range of —2.31 < [Fe/H] < —0.71 in 
the Victoria-Regina isochrones, as can be seen in Fig- 
ure [33] Although the data of Stream E were obtained 
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Fig. 25. — The Stream C and D regions shown zoomed-in (from Fig. [21]) with shces in [Fe/H] (-0.71 < [Fe/H] < +0.00 for left panel and 
-1.71 < [Fe/H] < -0.71 for right panel). 
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Fig. 26. — Contamination-uncorrected CMDs for two adjacent Suprime-Cam fields along the south-east minor axis, evalua ted with 
stringent criterion in SExtractor (see text). The superposed solid lines in the CMD show theoretical RGB tracks from VandenB erg et al.l 
(|200^) for an age of 12 Gyr, [a/Fe]= +0.3, and metallicities {left to right) of [Fe/H]= -2.31, -1.71, -1.14, -0.71, -0.30 and 0.00. 



under the worst conditions in our observing run, we can 
confirm well-populated metal-poor stars near the bright 
end of the M31 RGB. This population has a very similar 
CMD and surface brightness to those of the major-axis 
diffuse structure discussed in § 15.31 It also follows that 
Stream F shows a clear metal-poor RGB with a narrow 
metallicity range, as shown in the right panel of Fig. [33l 
This population has a very similar CMD to that of faint 
Stream A located at R 120 kpc along the south-east 
minor axis as found bv llbata et al.l (|2QQ7l ). 

The MDs corresponding to these two newly-found 
streams are shown in Figure [34l where we use the 



Victor ia-Regina theoretical isochrones for an age of 12 
Gyr and [a/Fe]= +0.3. In the north region of the M31 
halo, the number of Galactic dwarf stars exponentially 
increases as shown in Fig. [20l while the stellar density of 
the M31 h alo exponentially decreases with increasing ra- 
dius (e.g., Gilbert et al.' 2006). Therefore, because of the 
large uncertainty in correcting the disk contamination 
and the low statistics in extracting the intrinsic popula- 
tion of the M31 halo, we further impose the additional 
criterion of (V — /)o < 1.8 to construct the MDs, as was 
done for the MD of the major-axis structure (Fig. [3Q|) . 
The properties of the MDs clearly confirm the visual in- 
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Fig. 27. — MDs and cumulative MDs of the Stream C and D, 
using the Victoria-Regina theoretical isochrones for an age of 12 
Gyr and [a/Fe]^ +0.3 (as shown in Fig. I26|) . on the assumption 
that these substructures are as old as a diffuse, classical halo. 
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Fig. 28. — Luminosity functions of the Stream C and D around 
the RC. Solid line is the completeness-corrected LF with Poisson 
uncertainties, while the dotted line is the completeness-uncorrected 
LF. Red line is the best-fitted model with a simple model consisting 
of a Gaussian plus a straight line. We determine a Gaussian peak 
as the magnitude of the RC. Vertical dotted line denotes the 50% 
completeness limit. 
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Fig. 29. — Contamination-uncorrected CMD of the major-axis 
diffuse population, evaluated after stringent criterion in SExtractor 
(see text). 
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Fig. 30. — The MD and the cumulative MD of the major-axis 
diffuse structure, using the Victoria-Regina theoretical isochrones 
for an age of 12 Gyr and [Q;/Fe]= +0.3. We also select a bluer 
sample with (y — /)o < 1.8 to avoid the heavy contamination of 
disk dwarfs. 



spection of the above-discussed CMDs. Note that in this 
procedure we may miss some amount of the metal-rich 
stars with [Fe/H] > —0.7 because of this stringent crite- 
rion. 

6. THE GLOBAL HALO PROPERTIES ALONG THE MINOR 

AXIS 

6.1. Population Variation with Increasing Radius 

In Figures [35] and [36] we display the CMDs along the 
minor axis in order of projected distance from the M31 
center. After the background contaminations in the cat- 
alogs are subtracted by combination of DAOPHOT and 



SExtractor softwares in the same manner as our example 
above, we merge each catalog into final catalogs covered 
in 0.8-1.6 times field-of- views of a Suprime-Cam field and 
normalize them in the instrumental sky coverage with 
34' X 27'. It appears that the number density of the fore- 
ground Galactic disk dwarf stars predictably increases 
with decreasing galactic latitude, that is, in the northern 
direction from the M31 center and the number density 
of Galactic halo stars is practically constant along these 
galactic latitudes within the Poisson noise. This suggests 
that our method for the separation of galaxies and stars 
works successfully at the brighter magnitude. In con- 
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Fig. 31. — Luminosity functions of the major-axis diffuse struc- 
ture around the RC. 

trast, at the fainter magnitude of, for instance, /q ^ 23.5 
(/o > 23) for the SE 27kpc (41kpc) CMD, stellar objects 
seem to be over-subtracted because of the rapidly de- 
clining detection efficiency based on the above-mentioned 
stringent criterion. 

The two CMDs among the four south-east CMDs show 
the stream-like substructures (Stream C and D) as shown 
in § [521 The inner most CMD at SE 27 kpc has wide- 
spread RGB stars and it is very similar to that of the 
north-west counterpart field (NW 28 kpc), suggesting 
the existence of complex mix of metal-rich and metal- 
poor populations with —2.5 < [Fe/H] < 0.0 in the in- 
ner halo of M31. Both of the populations in the two 
symmetrically-located fields have the same mean metal- 
licity of [Fe/H]mean ~ -1.0. A KS test of the two cu- 
mulative MDs imposed by the additional criterion with 
{V — I)o < 1.8 (mag) in order to avoid the foreground 
contamination yields a probability of 12.3% that the two 
distributions are drawn from the same parent distribu- 
tion. 

In our Suprime-Cam survey of the M31 halo, we have 
obtained data of continuous fields along the north-west 
minor axis, facilitating the investigation of global halo 
structure. Furthermore, most regions in the north-west 
halo have yet been unexplored by any other telescopes. 
The continuous change of the CMDs along the north- 
west minor axis of Figure [36] suggests the existence of 
clear population gradients with distance from M31's cen- 
ter. The fraction of the metal-rich population with 
[Fe/H] ^ — 0.7 strongly decreases out to ~ 50 kpc in pro- 
jection, while the metal-poor stars with [Fe/H] < —0.7 
are definitely distributed in the outer part of ~ 50 kpc 
from the M31 center. This population variation identi- 
fied from the CMDs reflects a strong metallicity gradient 
with increasing radius in the M31 halo. In the more re- 
mote halo beyond ~ 50 kpc, the RGB sequence with any 
metallicity range is virtually unobservable, except in the 
two over dense regions such as Stream E and F fields. 

6.2. Surface Brightness Profiles 

Since iPritchet fc van den Berghl (|1994l ) derived the 
surface brightness (SB) profile of the inner dense parts 
of the M31 spheroid, many researchers have studied 
the extended, outer part s of the halo along the south- 
east minor axis (Durrell et al.l [2QQ1: llrwin ^nii [20051 : 
[Guhathakurta et al. 2005; Iba ta et "aL|[2007[ ). In contrast 
to the inner parts of the spheroid, it is possible to resolve 
the individual stars in its outer parts, because of their 
sparse density and of the sufficiently close distance of 
M31 to us. Thus the SB profiles in the spheroid's outer 



parts are available by directly counting individual stars. 
While previous studies have been limited to relatively 
bright RGB stars for the estimation of the SB in an ex- 
tended halo, we consider fainter stars, taking advantage 
of the light-gathering power of the Subaru telescope with 
an 8.2-m primary mirror. This suggests that the results 
obtained here using our deep Suprime-Cam data will be 
most robust against the effects of any contaminations 
over the fields. Furthermore, considering both the north- 
west and the south-east halo profiles simultaneously for 
the first time, we will arrive at more comprehensive halo 
structures than have been studied previously. 

We estimate the SB profiles using the RGB and AGB 
stars selected as (Vq^ Iq) < (24.8, 24.0), i.e., /-band 
magnitude brighter than the Red Clump (at Iq ~ 24.5) 
and colors redder than {V — I)o = 0.8, with not less than 
about 90% completeness limit. In addition we set Iq > 
20.3, i.e., fainter th an the tip of RGB stars of M31 (e.g., 
[Durrell et al.[[200Tl ), to exclude contaminations such as 
foreground st ars, when assuming the dis tance modulus 
of 24.43 from [Freedman fc Madorl (fl99Ql ). 

First, we divide a series of our surveyed fields shown 
in Fig. [1] into subregions with about 0.05 square de- 
grees. After extracting the secure RGB stars based on 
the above magnitude and color selection criteria, we di- 
vide them into the following two groups usin g the RGB 
model isochrones of [VandenBerg et al.l (|2006[ ) for an age 
of 12 Gyr and [a/Fe]= -hO.3: Metal-Poor group (MP) de- 
fined as -2.31 < [Fe/H] < -0.71 and Metal-Rich group 
(MR) defined as -0.71 < [Fe/H] < +0.00. This division 
is motivated by the form of the CMDs as obtained above 
(? l6.ip , showing a rather distinct feature below and above 
[Fe/H]~ —0.7 (corresponding to the metallicity of one of 
the Galactic clusters, 47 Tuc). 

Second, we convert the summed-up flux counts of se- 
lected stars to the SB in mag arcsec"^, and subtract the 
SB of the control field, for which foreground and back- 
ground contaminants are removed based on the statisti- 
cal method outlined in § I2.3.2[ The sources of noise in 
this method of estimating the SB arise from the Poisson 
statistics from the finite number of observed stars and 
subtracted contaminations. 

Figure [33 shows the reduced SB profiles of the MP 
(blue) and MR (red) groups. In this figure, two col- 
ored dotted lines having flat profiles present the sub- 
tracted background levels to obtain the reduced SB pro- 
files. In addition, we plot the SB in our observed re- 
gions of the GSS and the major-axis diffuse structures 
(green and grey) . The top panel of the figure denotes the 
SB profiles selected by our color cut of (V — /)o < 4.0 
used in the construction of the MDs in the previous sec- 
tions, whereas the bottom panel shows the SB profiles 
chosen by the bluer color cut of {V — I)o < 1.8 used 
in the construction of the MDs of the above-discussed 
outer streams. The latter color cut is motivated to avoid 
the heavy contamination of the disk dwarf stars. For 
guidance, a d ashed line corres ponds to the SB profile 
obtained by I bata et al.[ (|2007[ ), based on their survey 
for south outer parts of the halo; they identified stream- 
like substructures in some parts of our observed fields, 
whereby they derived the SB profile for the global halo 
by avoiding such disturbing substructures. This line is 
characterized by a power-law model, ^{R) oc with 
an exponent ofa = 1.91 ±0.12 kpc. We nominally fit 
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Fig. 32. — The Stream E and F region shown zoomed-in (from Fig. I21|) with shces in metaUicities (—0.71 < [Fe/H] < +0.00 for left panel 
and -1.71 < [Fe/H] < -0.71 for right panel). 
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Fig. 33. — The CMDs of Stream E and F, as derived in the same manner as Fig. 1261 The CMD of Stream E is made of combined catalog 
from 55 to 67 kpc, while the one of Stream F is constructed by summed-up catalogs from 92 to 105 kpc. The superposed solid lines in the 
CMD show theoretical RGB tracks from VandenBerg et al. (2006) for an age of 12 Gyr, [a/Fe]^ +0.3, and metaUicities {left to right) of 
[Fe/H]= -2.31, -1.71, -1.14, -0.71, -0.30 and 0.00^ 



this SB profile to our fields at ~ 30 kpc, which may 
lie outside the inner dense spheroid with a R^^^ law and 
trace an extended halo just beyond this component. It 
follows that this SB profile reproduces this inner part of 
an extended halo. At the brighter part of /iy < 33 mag 
arcsec"^, the decreasing trend of the estimated SBs in 
both the north-east and south-west parts is reasonably 
consistent with each other. We are also able to assess 
the faint stream-like substructures as identified in the 
population maps of Fig. [2T1 However, the absolute SB 
given in a dashed line is slightly underestimated in the 



inner region, which is dominated by metal-rich stars such 
as the W-shelf field. In contrast, at the fainter part of 
/iy > 33 mag arcsec"^, in particular, in the north-west 
region, the behavior of the SB profiles are somewhat dif- 
ferent. As discussed in § 16. H it is supposed that we do 
not detect both metal-rich and metal-poor stars in the 
outer region of > 80 kpc because we could not identify 
the metal-rich RGB with [Fe/H] > -0.7 on the CMDs 
(see Fig. [36]). Therefore, in the outer halo regions beyond 
R ^ 80 kpc we accept the SB profiles based on the more 
stringent color criterion of (V — /)o < 1.8, compared to 
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Fig. 34. — Metallicity distribution function and the cumulative 
distribution of the Stream E and F, statistically constructed in the 
same manner as the MD of the major-axis diffuse structure. 

those with {V — I)o < 4.0. Then, a prominent feature 
of interest is that the number of metal-rich stars with 
[Fe/H] > -0.7 in the north-west halo of M31 steeply 
decreases with increasing radius, compared to the south- 
east halo, because the south-east region within R ^ 70 
kpc may be polluted by substructures with metal-rich 
populations such as Stream C, D and the diffuse debris 
of the GSS as predicted by iFardal et all (l2007[ ). This 
picture is slightly inconsistent with our assumption that 
the region along the south-east minor axis at ~ 35 
kpc is free from any substructures, so we used the region 
for the control field through the analysis of the GSS in 
§ [3l However, we note that the effect of such substruc- 
tures does not significantly change the final result of our 
analysis for the GSS because the number of the contam- 
inated metal-rich stars at the R ^ 35 kpc field is only a 
small percent of those in our observed field of the GSS 
as shown in the right panel of Fig. [37l 

To investigate the spatial distribution of the M31 
halo in more detail, we attempt to fit the data on the 
radial profile (by weighting the Poisson errors in the 
star counts) of four theoretical models of the surface 
density profile. The power-law and exponential mod- 
els have been used to reproduce the surface bright- 
ness of some stellar systems (e.g., iBinnev fc Tremaind 
[1987). The NEW and Hernquist models represent- 
in^ the (dark matte r) density profile are derived from 
iNavarro et al.l (|1996f) (conv erted into the pr o jected for- 
mula by iBartelmannI 1 19961 ) and |Hernquist] () 19901 ) , re- 
spectivel y; this model choice is mo tivated by the simula- 
tions of Bullock fc JohnstonI ([2005). The dashed lines in 
Figure [38] show four different model fits to the selected 
minor axis regions, which are selected to probe a smooth, 
dynamically relaxed halo, where the contamination by 
substructures appears small. The inner halo, which may 
hold complex structures mixed by the debris of the GSS 
and/or the edge of the W-shelf structure, provides the 
upper limit to the radial SB profile of the minor-axis 
halo. 

Indeed, the outer halo profile appears remarkably flat 
in log-linear representation, thereby indicating essen- 
tially an exponential function. The black dashed line 
in Fig. [38] shows an exponential model fit to the data of 
MP group; we find a large exponential scale length of 
Rh = 18.8 ± 1.8 (22.4 ± 2.3) kpc for north-west (south- 
east) minor axis halo. We also show a projected Hern- 
quist model fit {thin short- dashed line) to these data; the 
best model gives a scale radius of Rh = 17.1 ±4.7 (31. 7± 
6.7) kpc, suggesting that the scale radius for north-west 



halo is comparable to that of lBullock fc JohnstonI (|2005l ) 
while that for the south-ea st halo is about a factor o f 
2 larger than pre dicted by iBullock fc JohnstonI (|2005l ) . 
Ilbata et all (|2007l ) also indicates the fiat density distri- 
bution as well as our south-east result, based on more 
outer and wider regions in the south quadrant halo out to 
140 kpc, because our south-east fields from 20 to 90 kpc 
are polluted by many substructures such as Stream B, C 
and D. 

These results, in particular for the north-west profile, 
are of great importance for understanding the nature of 
M31's halo, so these are to be assessed from a more care- 
ful analysis. We know that the number counts of back- 
ground galaxies rapidl y increases in the fainte r magni- 
tude range as shown in iKashikawa et aH (|2004l ). There- 
fore, there is concern about the misclassification of back- 
ground galaxies on the above calculated SB profiles. To 
test this effect, we re-calculate the SB profile along the 
north-west regions as before for the selected metallicity 
range, but this time we adopt an additional selection 
of stars: we further limit the halo sample with higher 
S/N, by imposing more stringent color-magnitude crite- 
ria of (Vb, Iq) < (24.3, 23.5). The results are shown 
in Figure [39] It is evident that these plots are qualita- 
tively and quantitatively identical to those based on the 
previous selection using deeper data. Also shown is the 
predicted behavior of the Galactic foreground contam- 
ination (with the same color-magnitude selection) as a 
dotted line in the figure. This profile of the contami- 
nation is nearly flat in this log-linear representation, so 
contamination cannot account for the observed profile. 
Thus, this smooth decline for the outer halo popula- 
tion of M31, whic h is comparable to the de nsity profile 
predicted by the B ullock fc Joh nstoiJ ('2005') model, is a 
robust result of this survey. We note that the density 
profile of the north-west minor axis halo derived here is 
steeper than that of the s outh-east minor axis halo de- 
rived bv Ilbata et al.l (|2007l ). This is because their sample 
is shallower than ours and they have underestimated the 
surface brightness in the denser, inner regions of the halo. 

6.3. Metallicity Distributions 

We derive the spatial variation of mean/median metal- 
licities along the minor-axis of M31's halo, by compar- 
ing the RGB stars on the CMP of each halo field to 
theoretical isochrones from "VandenBerg et al." (2006) in 
the same manner as previous constructions of the MDs. 
For adopted theoretical isochrones we assume the old 
and a-enhanced populations with an age of 12 Gyr and 
[a/Fe]= +0.3 like the Milky Way halo. Regarding our 
current method for the estimation of metallicities using 
these old iso chrones, it is noted that recent studies by 
iBrown et al.l ()2003l . I2006bl) suggest an intermediate age, 
i.e., younger than 10 Gyr, for about one third of the inner 
halo stars, in contrast to the Galact ic halo population . 
More recent studies with HST/ ACS ([Brown et al.ll2007l . 
l2008l ) have revealed the existence of truly ancient and 
metal-poor populations as well as younger populations in 
the inner halo of M31. While the age of these halo pop- 
ulations is yet to be fully understood, our method using 
the old theoretical isochrones provides appropriate lower 
limits on the metallicities of stars. However, as discussed 
in § 13.3. H for low-mass RGB stars with ages of 6-14 Gyr 
the metallicity derived from interpolation/extrapolation 
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Fig. 35. — The CMDs along the so uth-east minor axis of the M31 halo, 
tracks from IVandenBerg et al.l (|2006l ') for an age of 12 Gyr, [a/Fe 
-0.71, -0.30 and 0.00. 



The superposed solid lines in the CMD show theoretical RGB 
+0.3, and metallicities (left to right) of [Fe/H]= -2.31, -1.71, -1.14, 




Fig. 36. — - Th e CMDs along the north-west minor axis of the M31 halo. The superposed solid lines in the CMD show theoretical RGB 
tracks from IVandenBerg et al. (2006) for an age of 12 Gyr, [a/Fe]= +0.3, and metallicities (left to right) of [Fe/H]= -2.31, -1.71, -1.14, 
-0.71, -0.30 and 0.00. 



of bright RGB stars hardly depends on such an age dif- 
ference (Kali rai et al.|[2QQ6 b). 

We show their mean (median) metallicities as a func- 
tion of distance from the M31 center in the upper (lower) 
panel of Figure [JOl Several black filled circles present the 
average values of mean/median metallicities in our tar- 
geted fields, while numerous black open circles are those 
in the subfields, which are obtained by dividing a contin- 
uous field into the subfields with about 0.05 square de- 
grees. These subfields help us to examine a fine spatial 
variation inside a single Suprime-Cam field. Conversely, 
each filled black filled circle refiects the mean value of 
metallicities for the subfields. Vertical error bars of the 
upper (lower) panel show the associated standard devia- 
tion (interquartile interval) of the mean (median) metal- 
licities. 



These plots suggest that the stellar content of the 
M31 halo has no metallicity gradient within our observed 
range of fields. However, the spatial variation from 20 to 
50 kpc along the north-west minor axis appears to have a 
slight metallicity gradient. Our analysis of the CMDs dis- 
cussed in § [nH] robustly confirms the presence of the gra- 
dient. However, while the CMD analysis also predicts the 
lack of metal-rich stars in the more remote outer region, 
the figure shows the unexpected high fraction of metal- 
rich stars at > 80 kpc. This is because the derived 
mean/median metallicities in the outer halo are some- 
what overestimated due to the heavy contamination of 
the Galactic disk dwarf stars as discussed in § 16.21 There- 
fore, we re-calculate the mean/median metallicities by 
imposing a more stringent color cut of (V — /)o < 1.8, to 
ensure a minimal contamination from the Galactic disk. 
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Fig. 37. — The minor-axis profiles for RGB stars as derived from the MP and MR group, which are illustrated with blue and red 
points, respectively. Dotted lines denote the corresponding subtracted background levels. The error bars refiect a combination of Poisson 
statis tics and b ackground uncertainties. The dashed line shows a power-law profile with an exponent of 1.91 zb 0.12 for MP group derived 
by Ib ata et "all ((2007), for the sake of comparison. The distinction between the upper and lower panels refiects different color selection (see 
text). 



The results are shown in Figure [4T1 The mean/median 
metahicities in the south-east region show a systematic 
decrease of about 0.2 dex while keeping no spatial gradi- 
ent. In contrast, the variation in the north-west region 
based on the stringent color criterion shows now a strong 
metallicity gradient. 

We further test whether this strong metallicity gra- 
dient is a real feature or not, in addition to the above- 
confirmed morphological variation in the CMD along the 
north-west minor-axis. Figure [42] denotes the ratio of 



the signal-to-noise (S/N) ratio, plotted as a function of 
north-west minor- axis distance from the M31 center (in 
projection), where the S/N ratio means the ratio of the 
number of objects used to calculate the mean/median 
metallicity in each summed-up field to the number of 
the foreground sources in the scaled control field. The 
vertical axis of Fig. [42] means the ratio of the S/N ratio 
of (V — /)o < 1.8 selection sample to the signal-to- noise 
ratio of {V — I)o < 4.0 selection. Basically, because the 
vertical value is more than 1.0 in most regions of the 
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Fig. 38. — The north-west (south-east) minor-axis radial profile 
from the MP group is shown in upper (lower) panel. We reject 
the data containing the overdense regions as we have shown in the 
text. The thick long-dashed line is a power-law fit to these data 
with a = 2.17 ± 0.15 (1.75 ± 0.13), while the thick short-dashed 
line is an exponential profile with scale length of Rh = 18.8 zb 1.8 
(22.4 ± 2.3) kpc. The thin short-dashed line shows the best-fit 
Hernquist model, which has a scale length of = 17.1 ± 4.7 
(31.7 ± 6.7) kpc. In addition, with the thin long-dashed line, we 
show the NFW model for the halo mass profile with a scale length 
of Rh = 0.2 ± 3.5 (6.3 ± 3.9) kpc. 

north-west halo, the spatial variation of mean/median 
metallicities based on the stringent color criterion of 
{V — I)o < 1.8 is statistically more reliable than that 
based on a (V — /)o < 4.0 color cut. However, note 
that the absolute mean/median metallicities are clearly 
somewhat underestimated in the inner fields with a lot 
of metal-rich stars at < 40 kpc along the north-west 
minor axis. Eventually, we conclude that this test pro- 
vides the lower limit to the radial metallicity gradient 
photometrically detected in the north-west minor axis 
halo. 

7. DISCUSSION 

7.1. The Distribution of the Underlying Halo 

In this study, one of our major purposes is to clarify the 
general properties of the underlying, smooth halo com- 
ponent of M31. As mentioned in § (TJ recent large-scale 
photometric and spectroscopic surveys of the M31 halo 
have revealed the existence of the metal-poor and old 
populations (like a classical smooth halo in the Milky 
Way), which are buried in the i nner complex metal- 
rich and young s ubstructures (e.g.. iKalirai et aT]l2006bl : 
llbata et al.ll2007l ). By the term "smooth" , what we mean 
here is not necessarily that the component is perfectly 
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Fig. 39. — Same as the upper panel of Fig. (38] but for stars 
restricted in the small color-magnitude region Vq < 24.3 and 
Iq < 23.5, in order to more stringently remove the heavy back- 
ground contamination. Thus, these plots are only for the purpose 
of verifying the effects of such contamination on the derived surface 
brightness profile. The thick long-dashed line is a power-law fit to 
these data with a = 1.95zb0.18, while the thick short-dashed line is 
an exponential profile with scale length of Rh = 19.5zb2.3 kpc. The 
thin short-dashed line shows the best-fit Hernquist model, which 
has a scale length of Rh = 24.1 ± 7.0 kpc. In addition, with the 
thin long-dashed line, we show the NFW model halo mass profile 
with a scale length of Rh = 2.6 ± 4.5 kpc. 

spatially smooth, but instead that any substructures that 
may be present in the halo are below detectability with 
the current survey. 

The surface brightness profile as shown in § 16.21 is a 
good indicator to describe the radial extent of the halo. 
Since many studies of the density distributions of the 
Milky Way halo have already been made, it is useful to 
compare such results with our current studies of the M31 
halo. The Milky Way halo has been characterized by a 
steep power-law density profile p{r) oc r~^: its power-law 
index, 7, is in the range of 2.5 to 3.6. Some of the vari- 
ance in the value of 7 appear among different research 
groups, which can be attributed to the presence of sub- 
structures in the halo; previous studies have commonly 
been limited in sky coverage, so some local structures in 
the halo can affect a final result. Surveys covering signif- 
icant portions of the sky have been feasible o nly recently. 
For in stance, the most recent SDSS study bv lJuric et al.l 
(|2008f ). which made use of good distance estimates to 
halo stars near main-sequence turnoff covering 6500 deg^ 
of sky, found p{r) ex r~^-^ corresponding to a = 1.8 as 
a pow er-law inde x of the surface brightness profile. Sim- 
ilarly, iBell et dl\ ([2OO8 ) have found the stellar density 
distribution with p{r) cx r~^, based on about 4 x 10^ 
color-selected main-sequence turnoff stars contained in 
the SDSS DR5 catalog. Although in these SDSS surveys 
accurate investigations of the density profile are limited 
only within R ^ 20 kpc (or a t most 40 kpc, relied on the 
QUEST RR Lyrae sample of lVivas fc Ziiinl r2006) due to 
the uncertainty in the determination of the distance to 
remote stars, these results for the Milky Way halo are 
completely consistent with the present result of the M31 
halo, (X i^-2-i7±o.i5 (-^^^^ ^ _ 3 2)^ as derived for 

the north-west minor- axis halo out to R ^ 100 kpc. 

It is worthwhile to compare our current result of 
the radial profile of the M31 halo with other obser- 
vational studies dedicated to the outer halo. Our 
photometrically-estimated north-west profile is in rea- 
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Fig. 40. — Upper panel shows the mean metaUicity of north-west (left) and south-east (right) halo regions together with the associated 
standard deviation, plotted as a function of distance from the M31 center (in projection), while lower panel presents the median one together 
with the associated interquartile interval. Filled black circles present the average values of mean/median metallicity in the targeted fields, 
while black open circles denote those in each divided subfield (see text). Gray filled and open triangles present the mean/median metallicities 
of summed-up and individual fields of the major-axis diffuse structure, respectively. 
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Fig. 41. — Same as Fig. 1401 but for stars restricted to the small color region 0.9 < (V — I)o < 1.8, to ensure a minimal contamination 
from the Galactic disk. 



sonably agreement with the south-east profile of ex 
^-i.9i±o.i2 ]3g^gg(^ j^i^Q large south quadrant area out 
to 140 kpc (in projection) from the M31 center by 
the photometric CFHT/MegaCam survey of 'Ibata et al. 
(12007). In contrast, our resul t is slightly fl atter than 
iGuhathakurta et all ([2005f ) and lGilbert et all (|2006l) who 
have found an extended stellar halo in M31 follow- 
ing S(i?) (X R-'^-^ (or R-'^-^) out to - 165 kpc, 
based on the kinematically-selected RGB sample by the 
Keck/DEIMOS observation. However, we note that the 
former authors have not removed unresolved background 



galaxies which signiflcantly affect flelds under the worse 
seeing condition as discussed in § 12.3.21 and their survey 
is more than one magnitude shallower than our Suprime- 
Cam survey. The latter authors have estimated the radial 
proflle of the M31 halo derived from secure spectroscopic 
RGB sample in the remote outer halo, but it is based on 
only ~ 40 red giants to determine the outer proflle be- 
yond R = 60 kpc. 

Recent numerical simulations predict that most of 
the mass in the stellar halo originates in satellites that 
were accreted more than 8 Gyr ago (corresponding to 
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Fig. 42. — The comparison of the two S/N ratios, (S/N)i.8 [with 
a color cut of (V - I)o < 1.8] and (S/N)4.o [with (V - I)o < 
4.0], plotted as a function of north-west minor-axis distance from 
the M31 center (in projection), where the S/N ratio means the 
ratio of the number of the residual objects used to calculate the 
mean/median metallicity in each summed-up field to the number 
of the foreground sources in the scaled control field. 



at z :^ 1). The halo is smoothly built from the in- 
side out within the finite dynamical time scale, and be- 
comes the dominant contributor at 30-60 kpc and within 
(jBullock JohnstonI [20051 : lAbadi et an[2006h . Our de- 
rived surface density profile of the M31 halo in the ex- 
ternal view is in excellently agreement with the density 
profile of the accreted stellar halo; it is described as a 
Hernquist profile with Rh — 14 kpc within the radial 
range dominated by the an cient halo, as expected by the 
iBullock fc JohnstonI (|2005l ) model. In the meantime, our 
NFW fitting based on the observed SB profile returns 
several times higher concentration value (c = R^ir/Rh) 
than that derived from t he simulations based on dark 
matter distributions (e.g.. IBullock fc Johnstonll2005l ). as- 
suming the typical virial radius (e.g., i?vir = 200 kpc) of 
IVIilky Way- type galaxies. This fact suggests that the vis- 
ible baryonic matter like stars may well converge in the 
galactic center compared to the invisible dark matter be- 
cause of baryon dissipation during accretion processes. 

7.2. Implications for the Origin of Substructures 

Another major purpose of the current study is to probe 
the origin of various substructures as found in both the 
inner and outer parts of the 1VI31 halo. Their metallic- 
ity distributions as well as the spatial distributions of 
substructures are useful information to untangle the for- 
mation history of the halo because each substructure pre- 
serves such fossil records associated with star formation 
and merging/ accretion events. 

7.2.1. The South-East Minor Axis Substructures 

As shown in the previous subsection, the spatial dis- 
tribution of mean/median metallicity along the south- 
east minor axis of the 1VI31 halo remains approximately 
constant with increasing radius from the center. This 
is because our observing fields with Suprime-Cam coin- 
cidentally overlap some overdense stream-li ke substruc- 
tures, characterized by Stream D, C and B of llbata et al.l 
(|2007l ). which dominate the metallicity distribution of 
the stars in these fields. Stream D and C have sim- 
ilar high mean metallicity but clearly different CIVID 
morphology and metallicity distributions as shown in 



§ 15. 2[ suggesting that these two populations are differ- 
ent. Stream B brushes our outermost south-east field 
at ^ 90 kpc, as shown in the south quadrant map of 
llbata et al.l (|2007l ). This constancy in the mean metal- 
licity variation is consistent with the prediction by many 
previous studies, proposing that the halo is built-up 
through the dissipationless chaotic assembly fro m small 
subgalactic fragments (e.g.. iSearle fc Ziiinlll978f ). 

However, the abs ence of metahicity g radi ent confiicts 
with the result of Kali rai et al.l (2006b') and iKoch et al.l 
([2007), who found a strong metallicity gradient in the 
south-east minor axis discrete fields out to ~ 165 kpc, 
based on the kinematically-selected RGB sample with 
the Keck/DEIIVIOS spectrograph. The mean metallic- 
ities in the Subaru fields (SE2 a nd SE9) whi c h over - 
lap with the two fields studied in iKalirai et al.l (l2006bl ) 
(their a0[30kpc] and m6[87kpc]) are actually almost con- 
sistent with those given in their paper (([Fe/H]) ~ 1.1), 
although their work showed a metallicity gradient over 
larger radii out io R ^ 16 5 kpc. On the o ther hand, 
some discrepancies with the lKoch et al.l ()2007[ ) result are 
puzzli ng. These aut hors re-analyzed the same data as 
Kalirai et al.l (j2006bl ) but using a different method for 
metallicity determination and found a strong metallic- 
ity gradient, where the metallicity at ~ 90 kpc is 
Fe/H]= —2.4, i.e., more metal-poor than our result of 
Fe/H]~ —1. Although comprehensive understanding of 
this descrepancy is difficult only from our photometric 
studies. Stream B which distributes over our fields at 
~ 90 kpc may affect the determination of a global 
metallicity distribution, as mentioned above. To settle 
this issue, it is necessary to carry out a wider-field pho- 
tometric/spectroscopic survey over the large outer part 
of the M31 halo. 

In addition, we compare with the main-sequence fit- 
ting results based o n the HST/ ACS deep observations of 
iBrown et al.l (|2008l ) who found the trend that the mean 
age becomes simply older beyond 30 kpc in the south- 
east inner region. In fact, we have found that at 27 kpc, 
41 kpc (corresponding to Stream D), and 57 kpc (corre- 
sponding to Stream C), the mean ages are ~ 10, ~ 9.6, 
and ~ 9.3 Gyr, respectively, as derived from their MDs 
and RC niagnitudes using the t heoretically-computed re- 
lation oflReikuba et al.l (|2005h (see Fig. These re- 
sults suggest that there exist some inhomogeneities of the 
halo, which is yet on the way to the equilibrium struc- 
ture. 

7.2.2. The North-West Minor Axis Substructures 

The mean metallicity variation along the north-west 
minor axis, based on continuous Suprime-Cam pointings 
out to ~ 100 kpc, shows a strong spatial gradient with 
increasing radius. This is consistent with t he result of 
Kalirai et al. (2006b) and Koch et all (|2007r ). This spa- 
tial dependence of abundance distribution may confiict 
with the scenario of galactic halo formation via hier- 
archical merging, where successive merging events can 
erase any spatial gradients in metallicities as confirmed 
by numeri c al simu lations. For instance, the study of 
iFont et all (|2006bl ) predicts a gradient in [Fe/H] aver- 
aged over each full simulated halo of at most 0.5 dex 
over a few tens of kpc. 

It is yet difficult to clarify the physical origin of the 
detected metallicity gradient in the halo because there 
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exist few observations and simulations for the presence 
of a strong abundance gradient with increasing radius in 
the halo. A possible explanation for the current result is 
that the observed spatial distribution of metallicities may 
reflect just a snapshot of the dynamical evolution, which 
is on the way to the equilibrium structure characterized 
by a flatter metallicity gradient. In particular, the halo 
structure is dominated by transient substructures, which 
may largely affect the spatial distribution of metallicities 
instantaneously. Indeed, modern galaxy formation sim- 
ulations within the context of the standard ACDM cos- 
mology predict that the outer halo beyond R ^ 50 kpc 
is dominated by numerous faint arcs, shells and streams, 
as the relaxation time of these subsystems as found at 
large distances from the galaxy center is quite long at 
such location s, where even the dynamica l time is of order 
of Gyr (e.g., iBullock Jo hnston"2QQ5l: iJohnston etHI 
[2OO8). Furthermore, dynamical friction is at work, which 
acts more strongly on more massive subhalos, making 
them fall rapidly into the potential well, where they be- 
come disrupted and their contents are mixed into the 
evolving galaxy. We then observe a rather steep metal- 
licity gradient just after this event occurs. 

7.2.3. Dispersed Faint Substructure? 

Satellites with sufficiently long accretion times of more 
than several Gyr would contribute to the wide range 
of the stellar halo, so it is difficult to find any den- 
sity substructures in such well-mixed systems at the 
present day. It is noteworthy that previous studies found 
kinematically cold substructures based on the detailed 
Keck/DEIMOS spectroscopic observations in the south 
regions of the M31 halo, where no evidence for the pres- 
ence of substructures has been deduced from photomet- 
ric observations ^Ka lirai et al. 2006a; Gilbert et al. 2007; 
I Chapman et al. I [20081 ). These results suggest the possi- 
bility that there are also as-yet unrelaxed debris from 
disrupted dwarf satellites underneath the smooth halo 
at the north-west 30-40 kp c fields. Some A/'-body simu- 
lations (e.g., Bek ki fc Ch iba 2001) predicted that these 
optically faint substructures below our detected surface 
brightness levels may still exist in phase space, kinemat- 
ics. In this regard, we discuss about the possibility of 
the existence of these dispersed faint substructures in 
this subsection. 

As long as the surface brightness profiles (Fig. [37l and 
[38|) are concerned, there is no evidence for overdense sub- 
structures in the north-west and south-east regions at 30 
kpc. However, we find the age difference between both 
fields; the mean age of the stellar population in the north- 
west part (~ 7.5 ± 1.5 Gyr) is about 2.5 Gyr younger 
than that in the south-east counterpart field (~ 10 ± 1 
Gyr), as derived from their MDs and RC magnitud es us- 
i n^ th e theoretically-computed relation of Jlejkuba et 
()2005l ) (see Fig. [T6|) . In fact, we cannot insist that the 
two populations originated from the same stellar pop- 
ulation because the KS probability for this hypothesis 
based on their two cumulative MDs is rather low. There- 
fore, we expect that there might be a dispersed faint 
substructure in the seemingly smooth regions at = 
30-50 kpc in the north-west minor axis halo, which pre- 
serves the information of ancient accretion events. In 
addition, note that the age difference could be occurred 
by the difference in the distance of the same popula- 



tion substructures. To distinguish these much fainter 
substructures from the smooth inner halo we need addi- 
tional kinematic studies as well as some abundance in- 
formation (such as a-element abundance) as predicted in 
'Johnston et al.' (2008). 

Suppose that in the north-west region of M31 halo 
there exists such a well-mixed substructure with suffi- 
cient long accretion time (e.g., 9 Gyr), it may be possible 
to trace this accretion event from our photometric data if 
the event provides some specific stellar populations. If so, 
we expect a specific magnitude of the RC/AGBb feature, 
where a clustered distribution of stars in color-magnitude 
diagrams allows us to easily identify the presence of such 
stellar populations, compared to attempts to detect faint 
spatial substructures from already smeared out surface 
density distribution. Also we can estimate the mass of 
an accreted satellite as ~ 1 x lO^M©, which corresponds 
to [Fe/H] 

mean ^ — 1.0. Then, if wc suppose that this pro- 
genitor galaxy genuinely has an intermediate-age popu- 
lation, it is likely that such a galaxy possibly having in- 
terstellar gas has undergone a star formation driven by, 
for instance, the second or third collision or supernova 
wind about 7.5 Gyr ago. 

7.2.4. Comparison with Recent Numerical Simulations 

The recent models also predict that substructures that 
we can discover in the remote outer halo more likely have 
low surface brightness compared to the inner substruc- 
tures with high me tallicity and high surface brightness 
(e.g.. Bullock & Jo hnston 20051: IJohnston et al. 2008). 
In particular, Joh nston et al.l (|2008l ) predict that in a 
stellar halo of Milky Way type galaxies being built en- 
tirely from accretion within a ACDM context, the frac- 
tion of the system (of the order of 10% of its stars) are 
comprised of distinct substructures, i.e., a few to a dozen 
streams which are brighter than 35 mag arcsec"^ at pro- 
jected distances greater than 30 kpc from the host galaxy. 
Figure[43l shows the same cumul ative distribu tion derived 
from our M31 observations and llbata et al.l (|2007l) : the 
GSS, Stream A to F and the major-axis diffuse structure. 
The gray zone indicates plausible upper and lower limits, 
if simply extrapolated to the entire survey region. It is 
evident that both luminosity functions derived from the 
models and the observations are reasonably similar, sug- 
gesting that the outer halo of M31 as well as the Milky 
Way halo has been built up by recent accretion events of 
satellite galaxies. Furthermore, Figure SH plots the aver- 
age metallicity against SB for these substructures includ- 
ing our W-shelf data. Gray circles show the Stream A 
and B substructures of llbata et al.l (|2007l ). Since metal- 
licity uncertainties are just roughly plotted based on the 
error analysis in § 13.3. H in practice they may be slightly 
larger. Despite this, we find that higher SB tidal de- 
bris tend to be more metal-rich, which is confirnied by 
some recent numerica l simulations (e.g. JFont et al.ir2008l : 
IJohnston et al.|[2008l ). This result implies that higher SB 
debris may tend to come from the more luminous, metal- 
rich dwarf satellites and/or recent encounters. However, 
while more detailed theoretical approaches to the issue 
are required, it is also of great importance to explore 
more and fainter substructures in the larger region of 
the M31 halo using the next-generation wide- field im- 
ager such as Hyper Suprime-Cam, to set tighter limits 
on the halo formation history. 
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It is worth noting that recent h i gh-re solution A^-body 
simulation made bv lMori fc RichI (|2QQ8l ) provides an in- 
teresting scenario for the M31 halo: a merging event of a 
dwarf satellite into the halo of M31 gives rise to an elon- 
gated substructure similar to the GSS, the giant stellar 
stream found in the south-east part of M31's halo. Also, 
their model suggests that further stream-like substruc- 
tures emerge in the form of concentric shells around a 
galaxy after collisions with a disk component and sub- 
sequent relaxation process, and that a multiple large- 
scale shell system is finally formed in the outer region of 
the halo; the outermost large-scale shells have a radius 
of > 50 kpc and survive for at least 4 Gyr. If this is 
the case, then widely-separated shell structures would 
appear in the form of concentric shells around M31's 
halo. Although with this anticipation it is difficult to 
understand the existence of the clearly different radial 
metallicity gradients found along the north and south 
parts of M31, we actually have been able to identify two 
overdense substructures at the north-west counterpart 
fields. Considering their CMD morphology, it appears 
that the counterpart of Stream E is Stream B and that of 
Stream F is Stream A, although there is some positional 
shift of about 15-20 kpc. Since it is not yet clear that 
these counterpart streams indeed have the same popula- 
tions, we need to carry out further observations of the 
M31 halo with a wider and deeper imaging survey, in or- 
der to clarify the detailed stellar population and the true 
spatial extension of these substructures. 




36 34 32 30 28 

yu, (mag/arcsec^) 

Fig. 43. — The cumulative number of distinct substructures ob- 
served in the M31 halo, based on our work and Ibata et al. (2007). 
The gray zone indicates plausible upper and lower limits, if simply- 
extrapolated to the entire survey region. 



8. CONCLUSIONS 

We have investigated the structure and stellar popu- 
lations of the stellar halo in the Andromeda galaxy, us- 
ing deep and wide-field V and /-band images taken with 
Subaru/Suprime-Cam. From the analysis of the CMDs 
through comparing isochrones and globular cluster fidu- 
cials, the fundamental properties of each halo structure 
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Fig. 44. — The average metallicity again st SB for the subs truc- 
tures combining our M31 observations and llbata et al.l (|2007l ) (see 
text). Metallicity uncertaintie s (±0 .2 dex) are roughly plotted 
based on the error analysis in § 13.3.11 Higher SB tidal debris tend 
to be more metal-rich. 



identified here have been derived, such as distance, sur- 
face brightness profiles, metallicities and ages of the stel- 
lar populations. Our results for each of the GSS and 
other halo structures (substructures and smooth halo) 
are summarized as follows. 

8.1. The Andromeda Giant Stream 

We have identified an intermediate-age and metal-rich 
RC population as well as a broad RGB sequence in the 
Andromeda GSS. From the VI CMD and the LF, we 
have measured the distance to our observed stream field 
of (m - M)o = 24.73 ± 0.11 (883 ± 45 kpc) based on the 
/-magnitude of the TRGB and reddening E{B — V) = 
0.05. This result s uppor ts the previous prediction of 
iMcConnachie et al.1 ([20031 ) that their Field 5 (840 ± 20 
kpc) which falls near GSSl is at about 100 kpc further 
away along the line of sight relative to M31's main body. 

From the photometric comparison of more than 7000 
bright RGB stars in our stream field with the RGB tem- 
plates of Galactic globular clusters, we have obtained the 
metallicity ditributions (MDs) of the GSS for a statisti- 
cally significant sample. It shows the characteristics sum- 
marized as the peak metallicity of [Fe/H]peak ^ — 0.5, 
the average metallicity of [Fe/H]mean = —0.6 and the 
median metallicity of [Fe/H]med = —0.5. This re- 
sult is in excellent agreement with that of kinemat- 
i cally selected RGB stars in the same stream field 
(Guhathakurt a et al.| [2Q06). Furthermore, based on the 
stellar mass/luminosity- metallicity relation o bserved for 
the Local Group dwarf galaxi es shown in ICote et al.l 
(j2000r i and lDekel fc Wool (|2003f ). the mass of the plausi- 
ble progenitor is in the ran g:e of 10^ to 1 ^M(7^. which is 
consistent with the work of Mo ri fc Rich ' (2008 j. 

The feature near RC, showing an extended brighter 
part relative to its peak, is an unambiguous signature of 
an intermediate age population, which is also supported 
by the existence of an AGB bump. We have estimated 
the age of the stellar population in the GSS, by convolv- 
ing the MDs and the RC magnitudes with the theoret- 
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ical models from iRejkuba et aH (|2QQ5l ). As a result, we 
have found that the dominated population of the GSS 
has the mean age of ~ 8 Gyr, consistent with the main- 
sequence fitting results with the HST/ ACS observations 
bv lBrown et all i2006b). 

The GSS possibly originated from the accretion or mi- 
nor merger of a massive dwarf satellite within the mass 
range of 10^ to lO^M©, like M32, but previous studies 
suggests that M32 itself is probably not a progenitor of 
the GSS, as we already mentioned in Section 1. Consid- 
ering the finite dynamical time of such a mass accretion, 
the progenitor has to be bound to M31 several billion 
years ago as judged from the current residual substruc- 
ture. The first pericentric merging may well induce the 
third star forming event from the remaining gas, whereas 
30% of the GSS stars were provided in this event. 

8.2. The Andromeda Stellar Halo 

Our method for analyzing the CMDs of the GSS has 
also been applied to other large survey regions of the M31 
halo. Main results obtained from the outer halo survey 
are summarized as follows: 

• In the outer halo beyond R ^ 30 kpc, M31 con- 
tains a multitude of overdense substructures com- 
pared to the smooth radial profile, such as streams, 
arc and shells. We confirm the two substructures, 
namel y Stream C and D, discovered bv lb at a et al.. 
(|2QQ7l ). Furthermore, from our constructed stel- 
lar density maps we have found two as-yet un- 
known overdense substructures in the north-west 
halo, which we call Stream E and F in this work. 
It is yet unclear to what extent these substructures 
are spatially distributed. 

• We have derived the fundamental properties of 
these four stream-like structures in ad dition to the 
major-axis diffuse structure implied bv llbata et al.l 
(|2QQ7l ) and the Western shelf s een in the star count 
map of Ferguson et gj| (|2002l ). The derived prop- 
erties are summarized in Table [5l We have found 
that higher SB tidal debris tend to be more metal- 
rich, which is con firmed by some recent numerical 
simulations (e.g., iFont et al.l [20081 : iJohnston et al.l 
I2008D . This result implies that higher SB debris 
may tend to come from the more luminous, metal- 
rich dwarf satellites and/or recent encounters. 

• By separating regions showing obvious substruc- 
tures, we have found that the remaining area of 
the survey exhibits a smooth metal-poor stellar 
halo component in the north-west minor axis halo. 
This structure needs not to be perfectly spatially 
smooth, but the intrinsic inhomogeneities are be- 
low the sensitivity of this study. The smooth 
halo is vast, extending out to the radial limit of 
the survey, at = 100 kpc. The SB profile of 
this component can be modeled with a Hernquist 
profile as suggested by simulations, and the re- 
sulting scale radius of ~ 17 kpc is comparable 
to what modern halo formation simulations pre- 
dict. A power-law profile with ^{R) ex i^-2 i7±o.i5 
[p{r) (X 7--3-i7±o.i5j g^igQ fitted to the 
data. Our NFW fitting based on the observed 



SB profile returns a several times higher concen- 
tration value than that derived from the simu- 
lations based on dark ma tter distributions (e.g., 
iBullock Johnstoi] [2005h . assuming the typical 
virial radius (e.g., i?vir = 200 kpc) of Milky Way- 
type galaxies. This fact suggests that the visible 
baryonic matter may well converge into the galac- 
tic center compared to the invisible dark matter 
because of energy dissipation in the course of ac- 
cretion process. 

• We have detected metallicity inhomogeneities in 
different radial directions. The spatial distribution 
of mean metallicity along the south-east minor axis 
of the M31 halo remains approximately constant, 
while that along the north-west minor axis shows a 
strong spatial gradient with R out to 100 kpc. This 
fact suggests that the halo may be yet in the course 
of the dynamical evolution, leading subsequently to 
a fiatter metallicity gradient. 

In conclusion, the outer halo of the Andromeda galaxy 
as well as that of the Milky Way can be formed from dis- 
tinct stellar substructures deposited from over ~ 15 inde- 
pendent accretion events of subhalos with masses similar 
to typical dwarf galaxies, i.e., of order lO^-lO^M©. The 
Andromeda galaxy has quietly undergone the accretion 
of these intermediate -mass subhalos without the destruc- 
tion of the disk (e.g.. lKazantzidis et al.|[2Q08l ) since z ~ 1 
corresponding to ~ 8 Gyr ago. However, in the stellar 
halo of M31 there are probably much fainter substruc- 
tures below detectability with the current survey, as sug- 
gested by the existence of further relaxed substructures 
(e.g., slightly younger substructures discussed in § 17. 2p 
and kinematic substr uctures recently discovered in some 
spectroscopic s tudies (iKalir ai et al.|[2Q'06al : [Gilbert et al.l 
120071 : iChapinan et al. 2008). Therefore, we will need fur- 
ther observations with higher sensitivity to arrive at a 
better understanding of the formation of M31's stellar 
halo. 
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TABLE 5 

Fundamental properties of the substructures 



Name Distance /iy (typical) {[Fe/H]} (age) Progenitor mass"^ 



(Projected Distance from M31) 


kpc 


mag arcsec ^ 


dex 


Gyr 


Mq 


Giant Stream (~ 30 kpc) 


883 ± 45 


^ 29 


-0.74 


- 8 


- 109 


W-shelf ( < 25 kpc) 


798 ± 40 


- 30 


-0.75 


~ 8 


- 109 


Stream D (-^ 35-45 kpc) 




~ 31 


^ -1.0 


10 


- 109 


Stream C (- 50-60 kpc) 




~ 31 


^ -0.9 


~ 9 


- 109 


MA structure (~ 60 kpc) 




- 32 


- -1.2 


- 9? 


-108 


Stream E (~ 60 kpc) 




- 32 


- -1.3 




- 108 


Stream F (~ 100 kpc) 




- 33 


^ -1.7 


old? 


~ 10'^ 



* Based on the stellar mass-metallicity relation of lDekel Wool (|2003l ). 
the Ministry of Education, Culture, Sports, Science and Technology in Japan. 
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